VAPOR PHASE ACETYLATION OF ACETIC ACID 

TO VINYL ACETATE 


I. CHARACTERIZATION OF Zn(OAc)2 — 

II. KINETICS AND MECHANISM 

III. thermodynamic properties 


C SYSTEM 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

DOCTOR OF PHILOSOPHY 


by 

HARI BHAGWAN GOYAL 


to the 

DEPARTMENT OF CHEMICAL ENGINEERING 
INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 

MAY, 1976 



I j T. 

Jtoc. No. 


«) JL • 


2 3 SEP 1377 


CH£' 




(Vf^P 


i'Hi!. 




MY PARENTS 



CERa?IPICA[rE 


This is to certify that the work 'VAPOR PHASE 
ACETYLATIOF OP ACETIC ACID TO VINYI ACETATE, I. CHARACTERI- 
ZATION OP Zn(0Ac)2-C SYSTEM, II. KINETICS AND MECHANISM, 
III. IDEAL GAS THERMODYNAMIC PROPERTIES* has been carried 
out under our supervision and that this has not been 
submitted elsewhere for a degree. 


‘ 1 i L i 

A.P. ludchadker I 

Professor of Chemical 
Engineering 

Indian Institute of Technold 

Kanpur-208016, India ' 


( 35 ' 

R . D . Sr iva st ava 
Assistant Professor of 
Chemical Engineering 
Indian Institute of Technology 
Eanpur-208016 , Iiidia 


Date: May 


X 


1976 



[iii] 


ACKHOWLED&MMT S 

The author wishes to express his appreciation to 
Professor A,P, Kudchadker and Professor R.D, Srivastava 
for their inspiring guidance and advice, and for their 
close supervision and constant encouragement at all stages 
of this investigation. 

Dr. (Mrs) S,A.Kudchadker for her inspiring guidance 
in carrying out the thermodynamic studies. 

To Professor V., Srinivasan, Indian Institute of 
Technology, Madras for providing BET facilities. 

To Dr, L.M, Pandey and Mr, Raman of Defence Research 
Laboratories (Materials), Kanpur for their help in carrying 
out TGr analysis. 

To Professor D.KT, Saraf for his graciousness for 
providing constant inspiration during the work. 

To I4r. A,K. Agarwal, M.O. Garg, Rajeev Goel and other 
friends for their helpful cooperation from time to time. 

To the helpful cooperation offered by the staff of the 
Department of Chemical Engineering at Indian Institute of 
Technology, Kanpur particularly Mr. S. S. Chauhan and Ram 
Ashish Rai . 

To Mr, B.S, Pandey for his unfailing patience in typing 
the manuscript and Mr. D,S. Panesar for tracing the figures. 



his wife for her patience and encouragement 
throughout the course of this investigation. 

To his son for providing a welcome and refreshing 
change at home after a hard day’s work in the laboratory. 


Author 



CONTEHTS 

Page 

List of Figures . . . viii 

List of Tables ... x 

Nomenclature . . . xiii 

Synopsis ... xvi 

INTRODUCTION ■ ... 1 

SECTION I CHARACTERIZATION OF Zn(0Ac)2-C SYSTEM 

1.1 Introduction ... 4 

1.2 Experimental Methods ... 5 

1.2.1 Preparation of Catalysts 5 

(a) Preparation of Activated 

Carbon ... 5 

(b) Preparation of Catalj^-si® 7 

1.2.2 Physicochemical Characterization 8 

(a) X-ray Powder Analysis 8 

(b) Infrared ^ectroscopy 9 

(c) TG and DT Analysis 10 

(d) Surface Area, Pore Volume 

and Pore Size Distribution 10 

1.3 Results and Discussion 11 

1.3.1 X-ray Powder Spectra 11 

1.3.2 Infrared Spectra 13 

1.3.3 DT and TG Analysis 15 

1.3.4 Surface Area and Pore 

Distribution ... 18 



SECTION II 


SECTION III; 


[vi] 


Page 


: KINETICS AND MECHilNISM 


2.1 Introduction 

2.2 Experimental 

2.2.1 Overall Description 

2.2.2 Reactor 

2.2.3 Acetic Acid Vaporizer 

2.2.4 Plow Meters 

2.2.5 Condensers 

2.2.6 Purifier and Driers 

2.2.7 Experimental Procedure 

2.2.8 Analyst s 

2*3 Results and Discussion 

2.3.1 Catalytic Activity of Various 
Catalysts 

2.3.2 Effect of Zinc Acetate Concen- 
tration 

2.3.5 Comparison of Rate with. Surface 

iirea 

2.3.4 Effect of Acetylene-Acetic 
Acid Ratio in Peed 

2.3.5 Effect of ¥/P Ratios 

2.3.6 Effects of Diffusion 

2.3.7 Kinetic Models 

THERMODYNAMIC PROPERTIES 


24 

27 

27 

29 

52 

54 

54 

54 

55 

56 

58 

59 

59 

41 

45 

45 

47 

51 


5.1 Introduction 

for Tinyl 


64 

66 


5.5 


Gas Thermodynami 
of Vinyl Acetate 


c Properti-es 


5.3.1 Molecular Structure 

5.3.2 Vibrational Assignments 
5.3.5 Thermodynamic Puncti ons 


66 

66 

67 

72 



[vii] 


Page 


3.4 IdealcG-ae Pfe.ermodynamic' P'S'^perties 
of Cfotonaldehyde 

3.4.1 Molecular Structure 

3.4.2 Vibrational Assignments 

3.4.3 Thermodynamic Functions 

3.5 Acetone, Acetic Acid and Acetaldehyde 

3.6 Acetylene, Carbon Dioxide and Water 

3.7 Thermodynamic Properties of Mixtures 

3.8 Enthalpy, Gibbs Energy and 
Equilibrium Constants of Reactions 

3.8.1 Reaction Number 1 

3.8.2 Reaction Number 2 

3.8.3 Reaction Number 3 

3.8.4 Reaction Number 4 


74 

74 

75 

77 

80 

80 

81 

82 

82 

83 

83 

83 


CONCIUSION: AND RECOMMENDATIONS 

Conclusion 84 




Recommendations 

86 



REFERENCES 

87 

APPENDIX 

A; 

CH.lRilCTERIZilTION OF Zn(0Ac)2-C SYSTEM 

94 

APPENDIX 

B; 

KINETICS AND MECHilNISM 

100 

APPENDIX 

C; 

THERMODYNiMIC: PROPERTIES 

114 



. COMPUTER PROGRAM Cl 

131 



COMPUTER PROGRAM C2 

134 



+ + + •+*-{- 



+ + + 



[viii] 


LIST OF FIGURES 


Figure Page 

1 Differential Thermal Analysis Thermograms 16 

2 Analysis; Percentage Weight Loss of 
Zn(0Ac)2—0 ^stem as Function of Temperature 17 

3 Surface Areas as Function of Zn(OAc)p 

Composition ... 2o 

4 Relative Volume of Capillaries as Function 

of Capillary Radius ... p? 


5 Schematic Diagram of the Experimental Set-up 28 


6 

7 

8 

9 


10 


11 


12 


Reactor Details * ... 

Performance of Various Catalysts at 180°C 
and 25 per cent Metal Acetate by Weight 
on Activated Carbon ... ♦ 

Percentage Yield of Vinyl Acetate as 
Fxmction of Zinc Acetate Concentration at 
180^0; (¥/P) = 30; and (R)=3.0 

Effect of Reactant Mole Ratio, Acetylene 
to Acetic Acid, on Percentage Conversion of 
Acetic Acid at Various Temperatures, 
(¥/p)=30 ... 

Effect of Reactant Mole Ratio, Acetylene 
to Acetic Acid on Percentage Selectivity 
yinyl iicetate at Various Temperatures, 
(¥/E)=30 

Effect of Reciprocal Space Velocity on 
Percentage Conversion of Acetic Acid 
and Yield of Vinyl Acetate at 180°C and 
(H)=3.0 

Effect of Feed Velocity on Percentage 
Conversion of Acetic Acid at 180°C; 

(R)=3.0 and (¥/F)=50 ... 


30 


40 


42 

44 




50 



[ix] 


Figure 


Page 

13 

Acetylene Adsorption Model on Carbon 

56 

14 

Comparison of Observed' and Calculated 

Rates at 155°G and (W/P) = 30 

59 

15 

Comparison of Observed and Calculated 

Rates at 180® C and (W/P) = 30 

60 

16 

Comparison of Observed and Calculated 

Rates at 192^0 and (W/P) = 30 

61 

17 

Molecular Structure of Vinyl Acetate 

68 

18 

Molecular Structure of Crotonaldehyde 

S- Trans Porm ... 

76 

APPEEDIX 


A1 

P^re Size Distribution of Zn^^-000 

94 

A2 

Pore Size Distribution of Zn^'*'-005 

95 

A3 

Pore Size Distribution of Zn^'*'-010 

96 

A4 

Pore Size Distribution of Zn^'*'-015 

97 

A5 

Pore Size Distribution of Zn^''’-020 

98 

A6 

Pore Size Distribution of Zn^'^-025 

99 


+ + + + + 



-f + 4- 



LIST OF TABLES 


TABLE 


Page 

1 

’d' Values for Zinc Acetate 

12 

2 

Infrared Frequencies (cm~^) of Zn(0Ac)p-C 
Catalysts ... 

14 

3 

Surface Area, Pore Volume and ksr&CB.%'& Pore 
Radius of Various Catalysts 

19 

4 

Boiling Points of Organic Sabstances 

33 

5 

Retention Time and Area Mole Factors of 
Various Components . . . 

37 

6 

Role of Internal Diffusion 

48 

7 

Role of External Diffusion 

49 

8 

Rate Equations Derived Using Hougen-Watson 
Method , . , 

52 

9 

Comparison of Some In-plane Frequencies 

70 

10 

Adopted Fimdamental Frequencies of Vinyl 
Acetate (cm“l) ... 

71 

11 

Molecular Parameters of Vinyl Acetate 

73 

12 

Adopted Fundamental Frequencies of 
Crotonaldehyde ... 

78 

13 

Molecular Parameters of Crotonaldehyde 

79 

APPENDIX 

B 


B1 

Performance of Zn^'*’-025 Catalyst 

100 

B2 

Performance of Cd^'*'-025 Catalyst 

101 

B3 

2+ 

Performance of Ni -025 Catalyst 

102 

B4 

Effect of Zinc Acetate Concentrations 
on Per Cent Yield of Vinyl Acetate 

103 

B5 

Effect of R on Per cent Conversion of 

Acetic Acid and Per Cent Selectivity 
for Vinyl Acetate at 155®C 

104 





[xi] 

Table 


Page 

B6 

Effect of R on Per Cent Conversion of 
Acetic Acid and Per Cent Selectivity for 
Vinyl Acetate at 165°C ... 

105 

17 

Effect of Reactant Mole Ratio on per cent 
Conversion of Acetic Acid and Per Cent 
Selectivity for yinyl. Acetate at 180°C 

106 

B8 

Effect of Reactant Mole Ratio on Per Cent 
Conversion of Acetic Acid and Per Qent 
Selectivity for Vinyl Acetate at 1920C 

107 

B9 

Effect of Reactant Mole Ratio on Per Gent 
Conversion of Acetic Acid and Per Cent 
Selectivity of Vinyl Acetate at 199^0 

108 

BIO 

Effect of Reciprocal of Space Velocity 
on Per Cent Conversion of Acetic Acid 
and Per Cent Yield of Vinyl Acetate 
at 180°C 

109 

Bll 

Heat and Mass Transfer Effects 

110 

B12 

Comparison of Obssrved and Predicted 

R^tes at 155^^ 

111 

B13 

Comparison of Observed and Predicted 

Rates at 180®C 

112 

B14 

Comparison of Observed and Predicted 

Rates at 192°C ... 

115 

iiPPMDIX 

C 


Cl 

Vapor Pressure Boiling Point of Vinyl 
Acetate ... 

114 

C2 

Physical Properties of Vinyl Acetate 

115 

C3 

Martin-Hou Constants and the Departure 
Ettnctions for Vinyl Acetate 

117 

C4 

Ideal Gas Thermodynamic Properties of 

Vinyl Acetate ... 

119 

G5 

Ideal Gas Thermodynamic Properties of 
Crotonaldehyde ... 

120 





[xii] 

Table 


Page 

C6 

Ideal Gas ThermodyTaamic Properties of Acetone 

H 

CM 

H 

C7 

Ideal Gas Thermodynamic Properties of 

Acetic Acid ... 

122 

08 

Ideal Gas Thermod 3 mamic Properties of 
Acetaldehyde ... 

123 

09 

Ideal Gas Thermodynamic Properties of 
Acetylene ... 

124 

010 

Ideal Gas Thermodynamic Properties of 

Oarbon Dioxide ... 

125 

Oil 

Ideal Gas Theimiodynamic Properties of 

Wat er ... 

126 

012 

Enthalpy, Gibbs Energy, and Equilibrium 
Constant at Various Temperatures, Reaction 
Number 1 ... 

127 

013 

Enthalpy, Gibb H Energy,, and Equilibrium 
Constant at ^arieus Temperatures, Reaction 
Number 2 ... 

128 

014 

Enthalpy, Gibbs Energy, And Equilibrium 
Constant at Various Temperatures, Reaction 
Number 3 ... 

129 

015 

Enthalpy, Gibb'S Btiergy, and Equilibrium 
Constant at Various Temperatures, Reaction 
Number 4 ... 

130 


+ + -4- + -K 
-f + + 




[xiii] 



and b 


F 

& 

H 

P 

R 

R ■ 
S 
T 
V 
W 
X 
t 
d 

g 

It 


yrnMRRGLilTURl 


Mantin-Hou Constants 


Carboji 

Heat capacity 

Total flow rate (feed) 

Gibbs ^ energy 


ito.tlialpy 

ieat of reaction 
EGuilibri-um constant 
Total pressure 
Gas constant 





_ -1 « 4'’n 






^ ^ A 


Selectivity , entropy 
Temperature 


Volume 

Weight of the catalyst 
Conversion of acetic acid 


Yield 

Distance between^ each set 

the cirystal lattice 


of atomic planes of 


Gas 

Rate constant 




Adsorption constants 

^2 

Adsorption constant, rate constant 

m 

Order of diffraction 

P 

Partial pressure 

r 

Rate of reaction^ pore radius 

r 

Mean pore radius 

s 

Solid 

X 

Mole fraction 

Greek Symbols 

0 

Function 

V 

Pundamelital stretching frequency 

Q 

Angle of diffraction 

|il 

Microlitres 


Microns 


Wave length of x-ray beam 

5 

In-plane bend frequency 

Y 

Out-^f— plane bend frequency 

T 

Torsional frequency 

Subscrint 


A 

Acetic acid 

B 

Acetylene 

C 

Critical properties 

R 

Vinyl acetate 

S 

Surface 



0 


Property at 0°K 
a Assymetric 

s Symetric 

Super sc rj-pt 


o 


Property of ideal gas at 1 atm pressure 



[zvi] 


VAPOR PHASE ACETILATIOH OF ACETIC ACID TO VIHYL ACETATE 
I, Characterization of Zn(0Ac)2— C System, 

II, Kinetics and Mechanism, 

III. Thermodynamic Properties ' 
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FOR THE DEGREE OF 
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by 
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SYNOPSIS 

The thesis consists of three major sections. Section I 
deals with the physicochemical characterization of zinc acetate- 
activated carbon system. Carbon supported catalyst s-containing 
up to 25 weight per cent metal acetate were prepared by the 
impregnation technique. Section II deals with the kinetics 
and catalytic activity of zinc acetate-activated carbon system 
for the vapor phase acetylation of acetic acid to vinyl acetate. 
Calculations of the ideal gas thermodynamic properties of vinyl 
acetate and crotonaldehyde are dealt in section III. 

The vapor phase acetylation of acetic acid to vinyl 

% 

acetate over zinc, cadmium and nickel acetates supported on 
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activated carbon were examined at 180°C. Tbe catalytic 
activity increased as follows; Gd^'> Ui^'*' 

Supported zinc acetate on carbon catalysts containing 
up to 25 weight per cent zinc acetate have been studied by 
means of DTA, TGA, x— ray, IB. and BBI techniques in a wide 
temperature range (lOO® - 1000°C). Evidence is presented for 
a monolayer deposition up to about 15 per cent zinc acetate. 

The zinc acetate species, present in the dried catalyst did 
not interact with the carrier to form any stable carbon— zinc 
acetate complexes at the surface. The zinc acetate peaks 
were well defined. 

The heterogeneous interaction of acetylene and acetic 
acid with activated carbon based catalysts containing up to 
25 weight per cent zinc acetate have been investigated in an 
isothermal packed bed reactor at atmospheric pressure between 
155 - 200 C. The effect of several variables, molar feed 

ratio of acetylene to acetic acid, reaction temperature and the 
reciprocal of space velocity on the conversion and product 
distribution were examined by gas chromatography. 

Though 15 different mechanisms were postulated, the 
rate of reaction was most satisfactorily correlated by a 
mechanism of surface reaction between charged adsorbed acetic 
acid and acetylene, which assumes that the rate-controlling 
step was the irreversible charged adsorption of acetylene and 
acetic acid. 
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■where and Kg are adsorption equilihri'uin constants for 
acetic acid and acetylene, respectively* 

Complete ■vibrational assignments and calculations of 
the ideal gas thermodynamic properties of vinyl acetate and 
crotonaldehyde have been made for the first time from 298.15 
to 1500°K using the literature information on the fundamental 
frequencies and the molecular parameters. Internal rotational 
contributions due to CH^— group have been obtained for both the 
molecules using Pitzer and G-winn’s tables. 

For vinyl acetate the departure functions, that is, the 
difference between the real gas and the ideal gas property, 
have been computed using the Martin-Hou equation of state. 

As the pressure-volume-temperature data are not available for 
this compound, the Martin-Hou constants have been predicted from 
the critical constants and the vapor pressure-boiling point data 
of vinyl acetate. These departure functions in combination 
with the ideal gas properties will give the values of the real 
gas thermodynamic properties of vinyl acetate. It is hoped 
that in the absence of the experimental values, these predicted 
values will be of ■use to the design engineer. 


+ + + 



INTRODUCTIOIf 


Vinyl acetate an industrially impirtant monomer 
is polymerized into polyvin.yl 3,cetate which is used in 
latexes /• emulsion paints, adhesives, and various paper and 
textile coatings, or further processed into polyvinyl 
alcohol or polyvinyl butyral. Part of vinyl acetate produced . 
is also copolymerized with vinyl chloride, eth3?lene and other 
monomers . 

There are five different processes for making vinyl 
acetate (1-.4). These processes are differentiated primarily 
by their feed stock requirements, 

1. Acetylene-acetic aoid liquid phase process 

2. Acetylene-acetic acid vapor phase process 

3 . Acetaldehyde-acetic anhydride process 

4 . Ethylene-acetic acid liquid phase process 

5. Ethylene— acetic acid vapor phase process. 

Until. 1967, the acetylene vapor phase route accounted 
for most of the production of vinyl acetate. The ethylene 
route became operational in 1967 and because of the lower 
cost of ethylene as compared to acetylene offered attractive 
economics. However, the investment cost is stated to be 
about 50 per cent higher than for an acetylene-based plant. 
Electrical energy consumption is also higher. The picture 
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changed considerably in recent years. Also the coimtries 
who do not have adequate petroleum reserves but have large 
coal reserves need to seriously consider the acetylene 
route. 

In the vapor phase acetylene route, p-urified acetylene 
is mixed in excess with acetic acid in the vapor state and 
then fed to a fixed bed or a fluidized bed reactor using 
zinc acetate on activated carbon as the catalyst. The operating 
conditions are; temperature, 180-200°C; pressure, slightly 
above atmospheric. It is an exothermic reaction, the enthalpy 
of reaction being —23.4 kcal mol The selectivity of the 
zinc acetate catalyst is high and hence small amount of side 
products are formed, namely, acetone, acetaldehyde, crotonal- 
debyde, water and ethylidine diacetate. The separation of 
vinyl acetate from the reaction products is achieved either 
by fractional distillation or by liquid-liquid extraction. 

The yield of vinyl acetate varies from 92 to 95 per cent 
based on acetylene and from 97 to 99 per cent based on acetic 
acid (1-4), 

India imports all its requirements of vinyl acetate 
at present. Because of the availability of the raw materials, 
it was felt that the vapor phase acetylene-acetic acid 
technology should be developed in India. The present 
investigation is the first step ir. this direction. 
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The objective of the present investigation was to 
characterize zinc acetate-activated carbon catalysts, study 
the kinetics and mechanism Of the vapor phase acetylene - 
acetic acid reaction with the proposed catalyst, and to 
compute the required thermodj^-namic properties not available 
in the literature for the substances involved for process 
engineering calculations. Por the separation of vinyl acetate 
from other reaction products some vapor liquid equilibrium 
studies have already been carried out in this laboratorj'- (5). 
This information could then be used for modelling and scale- 
up which would be the second step in achieving the ultimate 
objective of providing the complete technical know-how for 
this process. 

The present investigation consists of three major 
sections. Section I deals with the physiccy^eijiieal characteri- 
zation of zinc acetate -activated carbon system. Carbon 
supported catalysts containing up to 25 weight per cent metal 
acetate were prepared by the impregnation technique . 

Section II deals with the kinetics and catalytic activity 
of zinc acetate-activated carbon ^stem for the vapor phase 
acet3'‘iation of acetic acid to vinyl acetate. Calculations of 
the real gas and the ideal gas thermodynamic properties of 
vinyl acetate and the ideal gas thermodynamic properties of 
crotonaldehyde are dealt in section III. 



SECTION I 


CHARACTERIZATION OE ZrL(OAc)2 -C STSTEM 
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1.1 INDEODUCTION 

f Catalysts based upon Zn( 0 Ac) 2 -C are remarkable for 
the acetylation of acetic acid to vinyl acetate. Fery little 
work has been reported in the literature concerning the 
physicochemical properties of these catalysts. All those 
reported have been carried out mainly by the Russian and 
Japanese workers (6-15). 

Surface area and pore structure of Zn(0Ac)2“C catalyst 
have been discussed by Yamada (9). Recently, Liverovskaya , 
et al (l5) have presented some data on the particle size and 
specific surface of various types of carbon, containing 
varying amounts Of zinc acetate. 

There appears to be no prior detailed study on the 
solid state properties of Zn(0Ac)2-C system. In order to 
contribute to the literature and to clear up some aspects 
concerning the nature of the active sites in relation to the 
catalytic activity, investigations have been carried out with 
a series of supported Zn( 0 Ac) 2 -C catalysts; proper choice of 
the support has allowed studies of compositions up to 25 
weight per cent zinc acetate. 

This investigation comprised the measurements of 
specific surface area, capillary structure, DT and TG analysis, 
infrared, and x-ray analysis. The acetylation of acetic acid 
has been chosen as the test reaction and its kinetics were 
studied (Section II ). 
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1.2 EKPERIMMTAL METHODS 
1.2.1 Pre-paration of (SJatalysts 

Por vinyl acetate manufacture from acetylene and 
acetic acid suitable catalysts include zinc or cadmium 
acetates, mercuric chloride, and others (2). These compounds 
may be supported on carriers such as activated carbon, wood 
charcoal, silica gel, activated alumina, etc. An improved 
technique for catalyst manufacture has been described by 
Bristol and Call (14). It comprised of dry mixing of zinc 
oxide powder and activated carbon, contacting the mixture 
with at least 4.5 mols of acetic acid and drying the catalyst 
by agitating under reduced pressure. They have claimed 
increased productivity and less linap formation. 

The overall procedure followed for mahing the catalyst 
was divided in two steps i 

(a) Preparation of activated carbon, 

(b) Preparation of catalyst, 

(a) Preparation of activated carbon ; 

The catalyst support activated carbon, was prepared 
in the laboratory using the steam activation technique. The 
basic idea of activating the carbon was to remove non carbo- 
naceous impurities from the surface and the pores, and thus 
providing more surface area and high adsorption power to the 


carbon 
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Coconut shells were used as the raw material. As 
desorihed hy HaSsler (15) it is a very suitable starting 
material for gas adsorbing granular carbon. The first step in 
activation was the 'cl^rring’ of shells. The term refers to 
the partial exidation by burning the shells, which is carried 
out in order to obtain a material with a high carbon content 
(about 90 per cent carbon) which may then be . subjected to 
activation and further oxidation 4lJ6^,* 

During the activation with steam the reaction between 
water vapor and carbon, even at temperatures as high as 700°G, 
involves first the adsorption of H 2 O on the carbon surface, 
followed by evolution of CO and H 2 . The CO is the most 
difficult energetically to remove from the surface. The 
overall reaction can be expressed as,«'’ 

H20(g) + C(s) ^^OOig) + H 2 (g) [ 1 ] 

and the rate expression f on reaction [l]> can he by 


Rate 


1+^2 ^jPb^O 



According to Strickland-Constable (17) the rate 
determining step in the mechanism is the slow desorption of 
CO from the surface. • 

In the laboratory the activation was carried out by 
suspending a bulk of carbon granules, packed in a 5 cm dia- 
meter and 25 cm long 316 stainless steel tube, in a furnace 
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heated up to 900°G. The steam was passed at the bottom of 
the stainless steel tube holder. The activation of each 
batch was carried out for about two hours and then the 
furnace was cooled and activated carbon was screened and 
stored. 

(b) Prenaration of catalyst : 

During the preparation of metal acetate on activated 
carbon catalysts, the equipments used were *r.otary Vacuum 
film evaporator and water ejector . For zinc acetate catalyst, 
calculated amount of zinc oxide based on the percentage of 
zinc acetate required in the final product, was dry mixed 
with 100 grams of granular activated carbon in the rotary 
vacuum film evaporator for one hour. ■ qlacial acetic acid, 
in weight ratio of two to one to zinc oxide, and water in 
approximate weight ratio of two and a half to one to acetic 
acid, were added. For example, to prepare 15 per cent Zh(0Ac)2 
sample the following proportions were required; 7*84 grams of 
zinc oxide, 15.5 grams of acetic acid, 100 grams of carbon 
and 40 grams of water. The whole mixture was agitated for one 
hour for completion of the reaction between acetic acid and 
zinc oxide to form zinc acetate on the surface of carbon. 

The whole mass was dried by agitating at 60° C under vacuxim. 

The vacuum was applied with the help of water ejector. The 
catalyst was further heated up to 150°G. 
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In the case of cadmiiim and nickel acetate catalysts, 
calculated amount of cadmium or nickel acetates based on 
percentage required, was dry mixed with 100 grams of activated 
carbon in rotary vacuum film evaporator. 0alculate4 amount 
of water was added, agitated, and the whole mass was dried 
by agitating under vacuum at 60°C followed by heating up to 
150°C. 

1.2,2 Physicochemical Characterization 
(a) X-ray powder analysis ; 

Each atom in a crystal has the power of scattering 
an x-ray beam incident on it. The sum of all the scattered 
waves in the crystal results in the x-ray beam being, in 
effect, diffracted from each allowed crystal plane. Every 
crystalline substance scatters the x-rays in its own unique 
diffraction pattern, producing a ’fingerprint’ of its atomic 
and molecular structure. The intensity of each reflection 
forms the basic information required in crystal structure 
analysis. One unique feature of x-ray diffraction is that 
components are identified as specific compounds. 

Virtually monochromatic radiation is obtained by 
reflection of x-rays from crystal planes. The relationship 
between the wave length of x-ray beam, the angle of diffraction, 
and the distance between each set of atomic planes of the 
crystal lattice d, is governed by the Bragg condition. 
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mA = 2d sin Q [5] 

where m represents the order of diffraction, Pirst order 
diffraction is, usually, used for structure study. 

The x-ray diffraction patterns were recorded using 
a General Electronics Diffractometer (XRD-VI generator, 

SPG-4 detector and SPG-2 diffractometer) with copper (E^) 
radiation. The diffractometer was operated with 2° diver- 
ging and receiving slits at a scan rate of 1* per minute 
and a continuous trace of intensity as a function of 2© 
was recorded. The accuracy of Bragg’s angle positions were 
checked using a permoquartz standard sample. Spectra of 
finely ground samples were run in the 2© range 10 - 80°, 

(h) Infrared snectroscp-py ; 

The range in electromagnetic spectrum extending from 
0,8 to 200 p is referred to as the infrared. Most organic 
and inorganic materials show absorption and, in all but a 
few cases, this absorption includes several characteristic 
wavelengths. In faet^ the infrared spectrum is one of the 
most characteristic properties of a compound. It provides 
a fingerprint for identification and a powerful tool for the 
study of molecular structure. 

Transmittance spectra for powder samples ( 4G0 mesh) 

were recorded on a Perkin-Elmer - 137 Grating IR spectro-. 
photometer using KBr disc. 
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(c) Itj and PI Analysis ; 

Thermogravimetric analysis (TG-A) involves changes 
in weight of a 'system -under investigation as the temperature 
is increased at a predetermined rate. Differential thermal 
analysis (DTA) consists of measuring changes in heat content ;, 
as a fiinction of the difference in temperature between the 
sample under investigation and a thermally inert reference 
compo-und, as the two materials are heated to elevated tempera- 
ture (or cooled to subnormal temperatures) at predetermined 
rate. In this manner enthalpic changes, such as melting, 
vaporization, crystallographic phase transition, or chemical 
changes, are detected from the end o-and exothermal bands and 
peaks that appear on the thermograms. The corresponding 
weight changes are detected by TGA. 

DTA measurements were carried out on a Dupont 900 
apparatus equipped with high temperature cell, using calibrated 
Platin'um/PlatinTim-13 per cent Rhodium thermocouples. TGA was 
performed on Gahn Thermo-Gravimetric apparatus, with programmed 
temperature increase, under atmospheric conditions, 

(d) Surface area, -pore volume and pore size distribution ; 

Surface area (SA) and pore structure {¥p and PSD) 
of catalysts are of major importance to catalysis since 
catalytic rates depend primarily on a-vailable active surface. 
Pore structure affects surface access, surface stability, and 
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resistance to poisoning, andselectivity as well as heat 
transfer. To correlate adsorption data and to determine 
the surface area, the multimolecular adsorption equation 

ef (BETT is offish Applied-, ^ 
Surface area of powdered catalysts were calculated 
from Ng adsorption isotherms using BET apparatus. Pore 
volume and pore size distribution were measured hy a mercury 
porosimeter, Carlo-Erha Model -70 for the 2000 atmosphere range. 

1.3 RESULTS AND DISCUSSION 

Reference to the various samples will be made by 
giving their composition as weight per cent Zn^t.025 

means a sample of zinc acetate supported on activated 
carbon having 25 per cent of the active element as Zn(0Ac)2* 

1.3.1 X-ray n’owder aoectra ; 

The prepared samples of zinc acetate were examined 
under x— ray to determine the purity. The 'd' values for zinc 
acetate are given in Table 1, calculated from x-ray powder 
diffraction patterns. The ’d’ values are in excellent agree- 
ment with ’d' values reported for zinc acetate (ASTM 1-0089). 
The x-ray diffraction pattern showed that the prepared 
samples have pure zinc acetate and there was no impixrity 
within the detection limit of the x-ray, - . 

The x-ray powder spectrum of the Z.n^-015,Zn^'^-020 
and Zj* —025 samples showed the presence of crystalline zinc 
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TABLE 1 ; »d* VALUES FOR ZINC AGBTATE 


A STM 1-0089 
'd* Values 

(S) 

Calculated 
’d' Values 

(A) 

7.50 

7.49S' 

4.55 

4.529- 

4.02 

3.92€; 

3.68 

3.693' 

3.46 

3.477 

3.28 

3.267' 

3.12 

3.145 

2.90 

2.949“ 

2.80 

2.735 

2.72 

2.622 

2.50 

2.495 

2.29 

2.16-5 
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acetate. The intensity of diffraction peaks were increased 
with the increase in the concentration of zinc acetate in 
the samples. In the 5-10 wt per cent samples the diffraction 
peaks of zinc acetate fell below the detection level. Powder 
spectra of all catalyst samples showed the absence of the 
characteristic diffraction bands of zinc oxide and hydrated 
molecules of zinc acetate, 

1.3.2 Infrared spectra 

Intensities of the various frequer.cies recorded by 
the transmittance techhiques are given in Table 2. Slight 
difference© are apparent in the ca.se ol pure compound, at 
variance with supported catalysts, which show distinct surface 
a.nd bulk structures. 

The stretching vibratiouval frequencies of the 0-0 and 
C-G bonds and the deformation vibration of the GO2 are listed 
in Table 2, which includes the corresppnding frequencies 
of the hydrated species. The spectrum of the prepared catalysts 
is very close to that of anhydrous zinc acetate (l8-20). 

X-ray spectrum agrees with this structure. 

In the highest zinc acetate concentration (25 weight 
per cent ZnCOAc)^) the transmittance spectra showed distinct 
absorptions at 1540, 1440, 950 and 694 cm~^, which are 
attributable to free zinc acetate present in the bulk. The 
decrease in the zinc acetate concenti*ation leads to the 



TABLE 2: INFRARED FREQUENCIES (Cm"-^) OF Zn(OAc)p-C CATALYSTS 
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absorptions which concentrate in a broad band in "the 1600 - 
400 cm~^ region. The supported catalyst of 5 weight per cent 
zinc acetate did not show any diagnostic vibrational frequency. 

1.3«3 DT and TG analysis 

A programmed heating at a rate of 15°C per minute 
between 25° and 1000°C was used. The DTA thermogram.s for the 
catalysts are shown in Figure 1, The thermograms revealed 
three regions of thermal activity^ (i) a weak endotherm peaking 
near 100°C very similar in appearance to the 100°C endotherm 
observed for the bare carbon, (ii) a strong and broad exotheim 
between 400° and 800°C peaking near 508°C, and (iii) a medium 
and broad endotherm (depending upon the composition of the 
catalyst) between 22cP and 360°C peaking near 280°C . The 
thermogram for pure zinc acetate exhibited a sharp endotherm 
of mediiun strength between 220° and 340°C peaking at 280°0 
and 320°C. The peaks at 280° and 310°C are due to the 
decomposition of zinc acetate to ZnO and Ac^O, and to ZnO and 
ac-otone, respectively. The results of the catalysts also 
agree with this; however, the intensity of the endotherms 
depend upon the composition of the catalyst. The peak in the 
500°C region represents burning of the carbon. 

TG- analysis was performed with programmed temperature 
increase at a rate of lO'^C per minute. The per cent weight 

t 

changes as a function of temperature are illustrated in Figure 2. 




DTA temperaturefC 




Fig. 1 '•Differential thermal analysis 







’ner ’T' ->3' cvimeJ 
Zr-ijAcyo-C sys 


ige-: weig'fli .:;ss 
ternp^ ■■ .: A.’ -- 




18 


The weight loss in the 200*^ - 320°C region is mainly hecause 
of the decomposition to zinc oxide. The weight loss above 
400*^0 is mainly due to burning of carbon* 

An analysis was also performed for pure zinc acetate. 
The thermogram was flat up to 180°0. Between 200°.- and 300°G 
a sharp weight loss was observed. The observed weight loss 
of about 60 per cent is in agreement with the calculated 
weight loss of 56 per cent (ZnO being the sole product). The 
results for the activated carbon showed no weight loss except 
above 400^0. 

The decomposition of zinc acetate has been studied by 
Edward and Hayward (19) by TGA and reported the formation of 
ZnO at 350°C as the sole product. The present results of 
DTA and TGA which are in complete agreement by themselves 
confirm this temperature range of decomposition, ZnO being 
completely formed below 350*^0, 

1.3.4 Surface area and pore distribution 

Results obtained by BET and mercury porosimetry are 
presented in this section. Surface areas (SA) and pore 
volumes (^ ) for all the catalyst samples are given in 
Table 3 together with the calculated values of the mean 
pore radii (r = 27^/SA). Values of SA for the samples are 
plotted in Figure 3 against values calculated on the basis 
of the contribution of activated carbon only. Figure 3 shows 



SURFACE AREA, PORE VOLUME AND AVERAGE PORE RADIUfj OT? 
■ VARIOUS CATALYSTS 
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variation of surface area with a reproducible discontinuity 
at 15 per cent and an almost linear decrease in the range 
of 0 - 10 per cent. 

The curves represented in Figure 4 clearly indicate 
that Zn( 0 Ac )2 deposited on the surface of carbon plugs at 
first small capillaries. As its amount increases, capillaries 
of greater and greater radii are closed and in the catalyst 
containing 15 to 20 per cent of zinc acetate, large capi- 
llaries dominate. 

It is necessary to point out that differentiation of 

curves represented in Figure 4 leads to the conclusion that 

greater values should be ascribed to capillaries whose 

o 

radius is below 100 A, This is in disagreement with the 

o 

value of the order of 10 A, calculated from the volume to 
surface ratio (Table 5). The latter value was however 
computed on the assumption that capillaries have tiCbe 
like shape . 

The pore size distribution curves for all the catalysts 
are given in Figures A1 - A6, Appendix A, 

An important point of discussion is whether zinc 
acetate is present on the surface of carbon as a monolayer, 
as a distinct phase in the form of clusters, or as interaction 
products with the carriers. Surface areas and pore volumes 
of the samples up to 15 weight per cent favor undoubtedly 
the monolayer deposition as they fit closely values calculated 




copaiory Todius. 
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on the basis of the contribution of support only (Figure 3). 
In catalysts with more than 15 per cent Zn(0Ac)2 increase 
in the mean pore radius and surface area suggests other 
processes. X-ray and IR failed to identify the formation 
of any well defined compound. It is possible that clustering 
as free ZnCOAc)^ may be occurring at these compositions. 


+ + + + -f 
+ 4 - + 
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II 


KINETICS AND MECHANISM 
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2»1 imaODUCTIOH 

During the last- 25 years, increasing attention has 
been given to the vapor phase catalytic acetylation of 
acetic acid to vinyl acetate. However, the entire source 
of literature on this subject is either in Russian or in 
Japanese. Among the various studies reported, there is 
wide disagreement in the reported data, and kinetics of the 
reaction are still not well understood (7,8,10,13,21-34). 

Purukawa (21) carried out kinetic measurements at 
200°C in a flow system an a zinc acetate-activated carbon 
catalyst. The order of the reaction with respect to acetylene 
was 1.65, and 0.65 with respect to acetic acid. From an 
investigation of the same reaction on a similar catalyst 
Janada and Vanko (26) obtained a second order equation. Also^ 
from a similar study at 120 mm and 160*^0 it has been suggested 
(23, 24) that the rate proceeds according to, 

^ ^CH,C00H 

r = 

(1+K Pqj£ QOOH^ 

3 

The reaction on the surface is the rate controlling step. 

Based on the theory that reaction of acetylene and 
acetic acid takes place via intermediate formation of 
complexes of metal salts i.e, catalysts with acetylene and 
acetic acid, Vasileva et al. (34) have examined the general 
scheme of the reaction mechanism by analogy with the hydro- 
chlorination of acetylene (39). ' 
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Zn(CH^COO )2 + CH^COOH ^;;^Zn(CH^a00)2‘CH^G00H 
ZnCGH^COO)^ + C 2 H 2 ^ Zn(GHjG00)2«G2H2 [6] 

Heaction between the complex of this salt with acetylene 
and acetic acid gives the following product ss 


Sn.(GH^G00)2*G2H2 + CH^GOOH — 5 ^ Zn(GH^GOO ) 2 + GH 2 = GHOOCGH^ 

[7] 

They assumed that the formation of a complex with acetylene 
by reaction [6] is the limiting stage of the process and 
derived a rate law given by, 


r = k. 




1 + 


\ ^GH^COOH 


[ 8 ] 


where k2 is the rate constant of the limiting stage, and 
is the equilibrium constant of reaction [5]. At temperatures 
above 180°G,. reaction is not retarded by acetic acid and 
equation [8] becomes a first order equations 


"■ = '"2 PCjHj [9] 

The chemisorption of acetylene is the rate determining 
step in the kinetics of the acetylene, acetic acid reaction 
(8,22,28,30,51). 

In addition to the heterogeneous mechanisms, a homo- 
geneous route involving free radicals is also believed to be 
involved ( 10 ) . 
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It is very probable that one of the causes of the 
discrepancy between the reported resiolts of kinetic measure- 
ments is inadequate isothermicity of the catalyst layer; in 
the given case this is an important factor, as the reaction 
is highly exothermic ( = -23.4 kcal mol“^) . These 

discrepancies may also be due to the mass transfer limitations, 
component concentrations, and high conversions. 

With increasing concentration of zinc acetate (30), 
the activity is proportional to the specific surface area and 
approximately to the cube of the zinc acetate concentration; 
but the surface area decreases exponentially with increase in 
zinc acetate concentration, and a decrease in activity at 
higher concentrations. The decrease of surface area with 
increase in zihc acetate is also reported by Mitsutani and 
Eominami (8). On the contrary Fakamura et al, (32) reported 
that the rate was independent of the amount of zinc acetate 
deposited on the catalyst. 

The body of results (8,13,30,32) from all these 
investigations has concurred in indicating singular;ities in 
catalytic activity as a f-unction of composition, thus justi- 
fying attempts to clear up relationship with its properties. 

In this work, the reaction between acetylene and acetic 
acid was studied in the temperature range of 150° - 200°C 
using zinc acetate supported on activated carbon catalyst. 

The size of the support, the metal acetate content of the 
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catalyst, and the reactor dimensions were made small 
enough to eliminate mass and heat transfer effects and 
permit a study of the intrinsic kinetics. In addition, 
an attempt has been made to correlate the results of this 
acetylation with the studies on the properties reported in 
Section I. 

2.2 EXPERIKEffTAl 

Preliminary runs were carried out in a 2.5 cm dia and 
12(0,, cm- , length packed bed reactor with zinc acetate supported 
on carbon catalyst. The conversion per pass obtained in 
the series of runs at 180°C was more than 50 per cent. With 
the success of these experiments it was then decided to 
study in detail the kinetics of the reaction employing the 
catalyst developed in this laboratory. The laboratory reactor 
was fabricated without giving any consideration to high 
conversion and yield. An overall description of the experi- 
mental apparatus, followed by detailed description of 
individual pieces is given below. 

2.2,1 Overall descrintion 

A general schematic diagram of the experimental 
equipment is shown, in Figure 5. 

The basic apparatus coiasisted of a fixed bed reactor, 
flowmeiers, control valves, purifier, driers, vaporizer, and 
condensers. All equipments were made of glass and interconnected 






KEY TO FiaURE 5 


1. Test Reactor 

2. Acetic Acid Vaporizer 

3. Product Condensers 

4. Soap Bubble Meters 

5. Silica Gel Driers 

6. Capillary Plow Meter 

7. Manometer 

8. Beedle Valve 

9. Acetylene Purifier 

10. Two ¥ay Stop Cock 

11. Acetylene Gas Cylinder 

12. Hitrogen Gas Cylinder 
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with, ground glass joints, polyethylene and silicon rubber 
tubings. The following capabilities and characteristics 
are provided fors 

(i) The flow rate of acetylene can be controlled 
relatively accurately for flows to less than 
1 per cent, 

(ii) The temperature of the reactor can be controlled 
to within + 0.5°G. 

(iii) The heat evolved in the exothermic reaction did 
not affect the temperature of the reactor. 

(iv) The output response from the thermocouples 
had .been measured accurately by means of a 
Honeywell potentiometer. 

(v) The product was completely condensed in the 
condensers. 

2.2.2 Reactor 

The jacketed reactor was 1.25 cm diameter and 60 cm 
long, provided with flexible preheating zone of maximum 
length of 16 cm. The preheating zone contained a 3l6 stain- 
less steel tubing packed with stainless steel foils to 
provide more heat transfer area. The stainless steel screens 
of 100 mesh size were attached at both ends of the preheaier. 
A diagram of the reactor is shown in Figure 6, 




Thermocouple 
well , 


Insufation 


Preheating 
zone ' , ' 


Fig. -6 - Reactor delaiis 
( dimensions in cms) 



32 


The temperature in the reactor was maintained by 
condensing vapors of various or^nic liquids. The substances 
used for various temperatures are listed in Table 4. 

The organic liquid was heated in a vaporizer. The 
vapors w'ere passed through the jacket, condensed in the 
condenser at the top, and refl^ized back. 

Heating by condensing vapors was found satisfactory 
in maintaining a constant temperature within + 0,5°C. 
Approximately two hours were required to reach the steady 
state temperature. 

The reactor was insulated by wrapping with approxi- 
mately 5 cms of glass wool and asbestos insulation to 
minimize heat losses to the surrounding. 

Temperatures at various points in the reactor including 
the catalyst bed and the heating liquid were measured with 
the help of iron-constantan thermocouples. The thermocouples 
were calibrated a.gainst a 25 ohm Leeds and Northrup platinum 
resistance thermometer and a resistance bridge assembly. 

The reactor was connected to the acetic acid vaporizer 
and the condenser (product) with the help of ground glass 
joints (Figure 5). 

2.2,3 Acetic acid vanorizer 

The vaporizer was a flask of one litre capacity 
and was heated by heating mantle. The desired temperature 
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BOILi p POINTS OP ORgAKIC -TOSTAFOEH 


temperature, ®C 

Substance 

155 

n-Hexanol 

165 

2-Methyl cyclohexanol 

ISO 

Aniline 

3-92 

Ethylene glycol 

199 

Y -Butyrolactone 



54 


of acetic acid vas obtained by supplying controlled voltage 
to the heating mantle. 

2.2.4 Flow meters 

The calibrated capillary flow meters (capillary 
size of 0,5 nim) and soap bubble meters were uOed to measure 
the flow rate of acetylene. 

2.2.5 Condensers 

Two double walled condensers were used in series to 
condense the product gases. Ice water was circuilat ed through 
the condensers to e.nsure complete condensation of the products 
This was confirmed by gas chromatographic analysis of the 
condenser exit gaa-es which contained mainly acetylene. 

2.2.6 Purifier and driers 

Acetylene was passed through a purifier and series 
of driers. The purity of the gas was checked by gas chromato- 
graphic analysis, 

and present in acetylene act as poisons to 
the catalyst. The purifier was used to eliminate these 
gases. 

The purifier contained a mixture of solutions of 
50 per cent CuCl^, 50 per cent PeCl^ and 55 per cent HCl in 
a volume ratio of 5i75l. This mixture could be reactivated 
by adding HCl to maintain 14 per cent HGl in the liquor. The 
gases from the purifier were dried with indicating blue 



35 


silica gel. The deactivated silica gel was replaced and 
reactivated hy heating in the oven at aro^md 100°G. The 
silica gel also absorbed acetone present in acetylene, 

2,2,7 Experimental nroced-ure 

The reactor was filled with fresh catalyst and the 
acetic acid vaporizer was charged with glacial acetic acid. 
The purifier and driers were filled with purifying solution 
and silica gel respectively. Nitrogen flow was started 
primarily to check the leak in the system and also to keep 
inert atmosphere while the system was heated during the 
experimentation. The reactor, heated by condensing vapors, 
and the acetic acid vaporizer were heated up to the required 
steady state temperature. 

Preliminary runs were taken in order to find the 
temperature of the acetic acid vaporizer at which it gave 
the required vaporization rate of acetic acid corresponding 
to a particular flow of acetyjlene. These runs indicated the 
temperature of the vaporizer to be maintained in a range 
of 100-117^0. 

The flow rate of acetylene was controlled with the 
help of a brass needle valve (Figure 5). ^mples were 
collected in sample bottles from the product condensers 
after every half an hour. The sampling bottles were kept 
in ice. In order to check whether vinyl acetate is 
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passing through the effluent , the uncondensed product 
gases were analyzed periodically, 

2,2,8 Analysis 

The analyses were performed by employing a Perkin 
Elmer 820 gas chromatograph provided with thermostated 
temperature control of detector, oven and injector, and 
temperature programming facilities. The output response 
of the chromatograph was recorded on a Honeywell recorder. 

The analysis of each sample reqtiired aroTund 25 minutes. 

A 300 cm long and 0,64 cm diameter 316 stainless 
steel tube, packed with 20 per cent FPAP on chromosorb P 
was used as- the analyzing column in the chromatograph. 

Nitrogen was used as a carrier gas (60 ml per minute flow 
rate). Temperature of the detector and the injector blocks 
were maintained at 250*^0 and 220°C, respectively. Temperature 
of the oven was kept at llO^G for acetylene, acetone, acetal- 
dehyde, -vinyl acetate, benzene, water, and crotonaldehyde 
elusion. For acetic acid the temperature programming was 
used from 110^ to 160*^C at a rate of 16°C per minute. 

Sample size of 2 pi was used. 

The retention time and area mole factors of' various 
components were found and are given in Table 5. 

During the runs, different feed ratios (acetylene/ 
acetic acid) were obtained by adjusting the acetylene flow 
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^BLE 5i I^TBETIOIT TIME AlH) AREA MOLE FACTORS OP 

VARIOUS OOMPOTJEWT.q 


Component 

Retention Time 
(Approx) , 
secs. 

Area Mole 
Factor 

Acetaldehyde 

90 

1.744 

Acetone 

120 

1.000 

Vinyl Acetate 

155 

1.669 

Benzene 

200 

— 

Crotonaldehyde 

290 

0.892 


Acetic Acid 8 minutes after the 

start of temperature 
programming 


1.604 
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rate in conjunction with, the temperature of the vaporizer, 
while the reciprocal of space velocity was varied hy changing 
the amount of catalyst and the total feed rate. 

2,3 RESULTS ilKD DISCUSSIOl 

Experimental data were obtained by means of an 

isothermal fixed bed reactor. The steady state was not 

only realized from the operating conditions but also from 

the product analysis. There was essentially no deactivation 

of the catalyst within 12 hrs. X-ray diffraction patterns 

showed that the used catalyst retained its activity. 

ilmong the byproducts, only acetone was in a measurable 

quantity. Acetaldehyde and crotonaldeh 5 rde were detectable 

in the gas chromatograph analysis, but the intensities of 

thermal conductivity were too small for accurate measurements. 

The effect of various variables, acetylene to acetic 

acid in feed (R), reaction temperature (T), the reciprocal 

of space velocity, (W/P)^on the conversion of acetic acid, 

the yield of vinyl acetate (VA), and the selectivity for 

vinyl acetate were investigated. The conversion (X), the 

yield of VA or acetone (Y), the selectivity (S), and the 

rate of the reaction (r) are defined as follows s 

conversion (per cent) 

- iQoles of acetic acid reacted ^ 

“ moles of acetic acid fed 
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yield (per cert) 

- n^ol^s of 7A or acetone -produced ^ 

“ moles of acetic acid fed ^ 

selectivity ( per cent) 

_ moles of Vil -produced 

moles of acetic acid reacted ^ 

rate of reaction 

_ moles of acetic acid reacted 
JhrTJ gm of catalyst) 

2.3.1 Catalytic activity of various catalysts 

file catalytic activities of zinc, cadmium and nickel 
acetates supported on activated carbon were examined^ The 
conditions were as follows: reaction temperature, 180°C; 
concentration of metal acetate by weight y 25 per cent, 
catalyst weight, 15 gm. The data are summarized in Tables 
B1 - B3, Appendix B, 

The yield of vinyl acetate against feed ratio (R), 

for various catalysts are shown in Figure 7. For the Zn 

sxib-groups the catalytic activity increases as follows: 

2 + 2 + 2 + 

Zfl ^Cd HI . Since zinc acetate on activated carbon 

gives the best yield, all experiments described below were 
made exclusively with zinc acetate. 

2.3.2 Effect of zinc acetate concentration 

The experimental data for the effect of zinc acetate 
concentration on yield are presented in Table B4, Appendix 3, 
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Figure 8 shows the per cent yield of vinyl acetate plotted 
as function of the concentration of zinc acetate. The 
reaction conditions were as follows; reaction temperature, 

180°C ; (W/F) = 30 gm-cat.hr mol^j (R) = 30. The yield 

passes through a maximum at the concentrs-tion of about 15 
per cent. This is in agreement with Rostovskii and Arbuzova(7) 
who have reported that the yield of vinyl acetate passed 
through a maximiira over a. catalyst concentration of about 
18 weight per cent ZnCOAc)^. 

Since 15 weight per cent zinc acetate on activated 
carbon ■sjas proven to give the maximum yield of vinyl acetate, 
all experiments described bel»w were performed exclusively 
with the catalyst of this composition. 

2.3.3 Gomoarison of rate with surface area 

It has been shown in Section I that the surface area 
decreases with the increase in zinc acetate composition up 
to about 10 per cent and then increases from 10 per cent to 
15 per cent . This compares well with the ’maximum rate observed 
at: 180°C at' a compositioir of 15 percent with the surface area, 
Rostovskii and Ushakov (30) have obtained a limiting 
value of activity which is proportional to the specific 
surface, and beyond this yalue the activity decreased abruptly. 
The maximum activity has been assigned to the unimolecular 
layer. 
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2.3.4 Effect of acetylene acetic-acid ratio in feed 

A study of the reaction rate of acetic acid ’^rer^s 
the feed ratio should give more precise information on 
the role of feed ratio In the conversion and selectivity 
of vinyl acetate. Reaction kinetics was examined in the 
temperature range 15^-200^0. The experimental data are 
presented in Tables B5 to B9, Appendix B, 

Figure 9 shows the effect of feed ratio (R) on the 
conversion of acetic acid over 15 per cent zinc acetate 
catalyst in the temperature range 155° - 200°G, for a (W/F) 
ratio of 30. The conversion of acetic acid increases with 
the feed ratio as well as with the reaction temperature. 

Figure 10 shows the effect of feed ratio (R) on the 
selectivity. It is seen that though the selectivity increased 
up to a temperature of 190°C at all (R) ratios, the change 
in the selectivity for the temperatures between 190°-200°C 
was very small. 

2.3.5 Effect of (W/F) ratios 

The effect of the reciprocal of space velocity (W/F) 
on the conversion of acetic acid, and the yield of VA and 
acetone were examined. The reaction conditions were as 
followsj reaction temperature, 180° G j (R) = 3.0. The data 
are summarized in Table BIO, Appendix B, 

The results of Table BIO are presented in Figure 11. 

It is seen that the per cent conversion and yield of VA 




f W/F ) j30 




*/o conversion IX] or yield (Y) 





'Q^ ^ conversion , ' 

• •% yield, vinyl', acetole 

• ' % yield ' aceto.ne' ' 


V '' 2 ' 0 ''/„ ,30 .40 50 ,:''''-',, 60 . V, , 70 ,'., 

,4,,^ ., ■' ' -Reoiproca! space vefacJly’ I W/F'). V' 

Fig. f1 -Effect of - reciprocaf space vefocjty on percentage 
conversion , of acetic acid and yieU vinyl acetate 
' iof 180“C and (Rti 3^^ 
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increased rapidly with the (¥/F) ratios studied, 

2,3.6 Effects of diffusiph 

The effects of diffusion were kept at a minimum by 
using high velocity of the gas through the catalyst. The 
reaction conditions were; (R) = 3.0, reaction temperature, 

- 180*^0 j (¥/F) = 30, The experimental data are shown in 
Table Bll, Appendix B, 

Figure 12 shows the effect of feed velocity on the 
conversion of acetic acid. The fair constancy of conversion 
obtained by changing the feed ratio while keeping (W/F) 
constant suggested that the diffusion of the gases was not 
controlling the rate. 

The role of internal diffusion has been checked by 
vaiying the size of the catalyst particle while maintaining 
the catalyst weight, (R), and (¥/F) constant. The internal 
diffusional resistance in the catalyst particles was practi- 
cally negligible as shown by fLnm gnif icant change in the 
reaction rate on varying the catalyst size between 2,8 to 1,4 
mm (Table 6) . 

In order to study the e ff ect of external diffusion 
it was necessary to increase the bed height while keeping (R) 
and (W/F) constant. The results are presented in Table 7. 

The rate was unaffected by doubling the bed height. Thus, 
the external diffusion in the catalyst bed was negligible. 
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TABLE 6 ; ROLE OF INTBRHAL DIFFUSION 

Effect of particle size of catalyst on rate of 
formation of vinyl acetate per hour per gram of 
catalyst at 180°C. 

¥ — 

j = 30.0 ; R = 3-.0 ; W = 15.0 


Run 

Fo. 


Average Particle 
Diameter, mm 


(gm mol Vinyl Ace tat 
formed hr~^ gm cat" 

X 10^ 


) 


ID-1 2.8 3.48 

ID-2 2.0 3.36 


ID-3 


1.4 


3.42 


H ffi 
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TABLE 7 : ROLE OP EITEMAL DIFFUSION 

Effect of doubling catalyst bed height on rate 
of formation of vinyl acetate per hour per gram 
of catalyst at 180^0. 

R = 3.0 


Run 

No. 

W, 

Weight of 
Catalyst , 
gm- 

F, 

Total Eeed 
mol hr~^ 

E 

gm cat hr 
mol”^ 

(gm mol Yinyl 
Acetate Formed 

hr“^gm cat“^)xl0' 

ED-l 

6 

0.6 

10 

5.20 

ED-2 


0.3 

20 

3.96 

ED-3 


0.2 

30 

3.70 

BD-4 

12 

1.2 

10 

5.20 

ED-5 


0.6 

20 

3.96 

ED- 6 


0.4 

30 

3.54 
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2,3.1 Kinetic models 

Several kinetic models can be postulated based on 
the chemisorption behavior of zinc acetate— C system. The 
mechanism of the solid catalyzed gas reaction, in general, 
may be the mass transfer of the reactants, or products, or 
adsorption of the reactants or desorption of products, or 
chemical reaction between adsorbed molecules at the catalyst 
surface. 

Assuming that the various acetylation processes 
involved have nearlj'- the same functional dependence on p« , 
^CH COOH* expression for the rate takes the forms 

Rate = k 0 -CE^COOE* %A^ 

Based on the conventional langmuir-Hinshelwood (35) 
mechanism a kinetic analysis of the experimental data was 
made using the approach suggested by Hougen and Watson (36). 
In this method, various mechanisms, which might control the 
reaction are postulated, based on a single site (Table 8a) 
and two types of active sites (Table 8b), 

Neither zinc acetate nor activated carbon show any 

appreciable catalytic activity by themselves, nor do they 

chemisorb acetylene. The heat of adsorption of acetylene 

on active carbon is only 4,2 kcal mol“^. It may be explained 

as followss since activated carbon is mainly composed of 

0 

graphite, a lattice spacing of 1.42 A could be chosen as the 


i ■ ,i. 



TABLE 8 ; RATE EQUATIONS DERIVED USINO HOUGEH-WATSON METHOD 
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characteristic distance between carbon atoms in the graphite 
structure (37). If acetylene adsorbs by breaking a bond 
followed by bond formation with the carbon atoms of graphite 

(ethylene structure), then the optimum carbon carbon distance 

o , . 

is 2.88 A for unstrained adsorption (38). Since the 

o 

distance of 2.84A exists in the graphite lattice, it may 
be possible that acetylene adsorbs on graphite by acquiring 
an ethylenjc structure. However^ the mode of adsorption and 
the energies involved would be expected to preclude easy 
mobility. A schematic model of acetylene adsorption on 
carbon is shown in Figure 13. 

It is generally accepted that carbon acts as a Lewis 
acid, defined as an electron acceptor, and the zinc acetate 
as Lewis base, as the electron donor (l). The transfer of 
electrons from donor to acceptor would result in ionized 
species. It was assumed that this surface consists of two 
types of active sites. On these active sites acetylene and 
acetic acid are selectively adsorbed. 

The heterogeneous catalytic reactions traditionally 
described as, 

(a) diffusion (external and internal) of 
reactants and products, 

(b) adsorption and desorption of reactants 
and products (physical adsorption and 
weak chemisorption), 
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and (c) surface reaction 

should be modified to take place as follows: 

(i) electron transfer from or to the reactants, 

(ii) homogeneous reaction (formation of activated 
complex and rearrangement of the ionized 
species on the surface), 

and (iii) electron transfer from or to the products. 

Writing the stoichiometric equation for the reaction 

as, 

C2H2(g) + CH^GOOH(g) ^H2=CH00CDH^ (g) [H] 

following the usual procedure of derivations (42) , and 
neglecting the term Py^/K, since E, equilibrium constant for 
the vapor phase aqetylation of acetic acid was very large. 
(in E > 12 ; temperature 150°-200°C), equations for various 
postulated mechanisms, relating the rate of reaction were 
derivedf (Table 8a and 8b), and tested for best fit of the 
data. 

Linear regression analysis of Hougen and Watson 
models to obtain the least squares fit values of the kinetic 
constants (43). were performed on IBM 7044 system. The 
program prepared used centralization and correlation matrix 
ifor minimizing the round off errors. The linear regressions 
of the rate equations were made at each temperature, Ai i 
possible combinations of the various adsorption terms were. 
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examined in order to allow for various negligible adsorption 
terms. Rate expressions were eliminated from further consi- 
derations when any of the converged adsorption constants 
became negative. 

The data were most satisfactorily correlated by 
mechanism 14 (Table 8). The square of the multiple 
correlation coefficient (R^) was found to be above 95 per 
cent for this model at all temperatures. The mechanism is 
based on the surface reaction between ionized acetic acid 
and acetylene. 


rate 


K 


■A Pa 


. X /PAx 
1 + ( IT 

Pr 


[ 12 ] 


where and Eg are temperature dependent and can be expressed 
by the following equations. 


log E^ = 6.305 - [13] 

log Eg = 0.33 - [14] 

These eonstante also fit the rate data. Table B12 - 
B14, Appendix B. The plots of the experimental rate versus 
that predicted by the above expression are shown in Figures 14, 
15 and 16 for the temperatures 155^, 180° and 192°C. The 
maximum deviation between the predicted values and the 
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experimental ones have been found to be + 7 per cent, 
except at temperatures above 190*^0. 

The model presented in this work differs from that 
given by Yasileva, et al, (34). In their study, the rate 
of reaction on zinc acetate on carbon was first order with 
respect to acetj^iene and independent of the partial pressures 
of acetic acid or vinyl acetate above 180 °q^ In analogy with 
the hydrochlorination of acetylene (39) a rate expression 
was presented. It has been well established that neither 
zinc acetate nor activated carbon show any catalytic activity 
by themselves, even at relatively high temperatures (180° - 
218°C). This distinguishes it very clearly from the hydro- 
chlorination of acetylene, which is catalyzed by one of the 
components. Careful examination of their results was not 
possible since most of the original data were not presented. 
Prom the discussion of various mechanisms it is 
obvious that linear regression applied only to rate data will 
not be sufficient to determine a definite mechanism. The 
assumptions which are made in the development of the particular 
rate expression must be experimentally verified. 

It is interesting to note that the form of rate 
expressions derived in Table 8 is similar to those derived 
by Mann and coworkers (40,41) for the oxidation of olefins 
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over bismuth molybdate where the surface was considered 
to consist of two types of active sites on which olefin 
and oxygen are selectively adsorbed. 


+ + + + + 


+ + + 
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THBRMODnTAMIO PROPERTIES 
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3.1 INTRODUCTION 

A detailed process engineering study of vinyl acetate 
manufacture requires the kinetic rate equation, vapor pressure- 
boiling point data, thermodynamic properties of the real gas 
at the process conditions, and the vapor- liquid equilibria 
of the multicomponent mixture comprising the reaction proudcts. 
In Section II, the rate expression has been reported based 
upon the results of the present investigation. 

The vapor pressure-boiling point data for vinyl acetate 
are reported by Capkova and Pried (44) and are given in 
Table Cl, Appendix C along with the Antoine constants and 
the enthalpy of vaporization. 

The real gas thermodynamic properties and the experi- 
mental pressure, volume, temperature data for vinyl acetate 
are not available in the literature. Hence we had to take 
recourse to some reliable predictive methods for obtaining 
an equation of state for vinyl acetate. The Hartin-Hou 
equation of state (45) was selected for this purpose and the 
Martin-Hou constants were computed according to the procedures 
outlined by Martin and coworkers (45, 46, 47). The departure 
f^inctions that is, the difference between the real gas 
property and the ideal gas property, were then derived using 
the following expressions involving the Martin-Hou constants. 
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Enthalpy content 

' P P 

= J dP = j [V-T{^)^(-^)’ ]dP 

o o 

[15] 

Entropy content (S-S°)^ 

P 

(S-S°)^ = -^J [(-||)y(-||);^]p + R In ? 

° [ 16 ] 

Gribhs energy content (G— 0-°)^ 

_ , 

(a-G°)^ = - J [ ^ - Y] dP + HP In ? [17] 

o 

Phe thermodynamic properties of the real gas is the 
Slim of the quantity for the ideal gas and the departure 
function calculated above. This section deals with the 
departure functions for vinyl acetate computed by the procedures 
outlined above. The ideal gas thermodynamic properties for 
vinyl acetate and crotonaldehyde are not available in the 
literature. Hence it was necessary to calculate these 
properties using the literature spectroscopic information and 
are reported in this section. 

The other properties for vinyl acetate are reported 
in Table 02, Appendix 0 (l,3) 
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3.2 m DEPARTURE gUWGTIOFS FOR YISTYL ACETATE 

The departure functions were calculated using the 
Martin— Hou equation of state given below; 


log 


■p 

~Y--b 



A. 


B.T 

i 


+ 


C.e^^/^c 




[18] 


Here to A^, to B^, to C^, k, b are the Martin-Hou 
constants and could be determined from the critical constants 
^c* ^c^ "^c’ point on the vapor pressure curve. 

Using these properties for vinyl acetate as reported in 
Table C2, Appendix C, the Martin-Hou constants were computed. 
Bor this purpose a computer program was written according to 
the procedures given by Martin and coworkers (45,46,47). 
Equations [15 to 18] were used to compute the departure 
functions at selected temperatures and voliimes in the vapor 
region and these along with the Martin-Hou constants are 
reported in Table 03, Appendix C, The computer program used 
for this purpose is reported in Appendix C. 


3.3 . IDEAL GAS THERMOBYIIAMIG .PROPBFMIES OE YIMYL ACETATE 
3.3.1 Molecular structure 

Two planar conformations are possible for vinyl acetate 



Gis 


Trans 
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where R and R’ are methyl and vinyl groups, respectively. 

As no microwave studies have been carried out for vinyl 
acetate, the molecular conformation is not known. However, 
from infrared studies (48,49) it is found that no rotational 
isomers are present, Rao and Curl (50) have established cis 
confonnation for vinyl formate from their microwave spectral 
studies. Crowder (51) interpreted the infrared and Raman 
data for vinyl trifluoroacetate in terms of one conformer of 
Cg symmetry, The review article by Jones and Owen (52) makes 
it very clear that vinyl acetate exists as cis conformer with 
Cg symmetry as shorn in Figure 17. Here the carbon oxygen 
skeleton has a planar structure and one hydrogen of the methyl 
group is in plane while the other two are symmetrically 
arranged with respect to it. The values of the structural 
parameters adopted in this work are given in Figure 17. These 
are taken from the parameters of similar molecules, namely 
vinyl formate (48), divinylether (53), methyl acetate (53,54) 
and acetic acid (55). 

3.3.2 Vibrational assignments 

Kotorlenko et al, (49) reported assignments for this 
molecule only from infrared measurements. Feairheller and 
Raton (48) thoroughly studied the infrared and Raman spectra 
and reported 19 in-plane and 9 out-of-plane fundamental 
frequencies. The missing in-plane frequency is and 
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the one out of plane is the torsional frequency due to 
methyl group, 

McManis (56) studied the infrared spectra of series 
of vinyl esters and assigned 464 cm“^ to the 6„ « „ mode. 
Crowder (5l) reported a similar value (465 cm“^) for this 
mode in the case of vinyl trifluoroacetate. Assignments of 
the three lowest in-plane frequencies of similar molecules 
are given in Table 9, It was felt that the two lowest in— plane 
frequencies reported by Peairheller and Eaton (48) should be" 
assigned. Therefore 462 cm~^ to the mode and 406 cm“^ 

to the mode have been assigned. This means that 

*^G-0-C been observed by Feairheller and Eaton, The 

missing 6 q__q__q frequency is found to be 306-523 cm”^ in 
aliphatic esters (57), 303 cm"*^ in methyl acetate (58) and 
295 cm”^ in vinyl trifluoroacetate (51). The average value 

299 cm has been- ad-opted in the calculations. 

The out-of -plane missing torsional frequency adopted 
in this work is 100 cm ^ which is transferred from methyl 
acetate (59). In acetic acid this value is reported to be 
93 cm“^ by Haurie and Novak (60) and 101 cm“^ by Fukushima 
and Zwolinski (6l). The fundamental frequencies adopted in 
this work are given in Table 10, 
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TABLE 9 s COMP ARI SPIT OF SOME lE-PIAME PREQUEMCIES 


Mode 

Methyl Acetate 
(58) 

Vinyl Trifluoro— 
acetate (5l) 

Vinyl Acetate 

Ref.(4S) 

' Work 

^ 0 = 0-0 

- 

464 

- 

462 

^ 0 - 0=0 

429 

356^ 

462 

406 

^0-0-C 

303 

295 

406 

299 

Low iralue because of 

mass effect of CP^ 

3 

group. 
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table 10; ADOPTED EOUDAMEFTAL EREQUBHCIES OF VINYL ACETATE 


In-plane frequencies (a*) 

1 

3120 

'iConp 

11 

1353 

5(CH) 

2 

3090 

jkcH) 

12 

1291 

^(C-C) 

3 

3040 

VsCCHj) 

13 

1217 

U(coc) 

4 

299 € 

\(0H3) 

14 

1018 

5(GH2)rock 

5 

2940 

■VeCOHj) 

15 

972 

6(CH^)rock 

6 

1760 

V(c=o) 

16 

847 

"j^^COC) 

7 

1644 

'))(o=o) 

17 

637- 

6(002) 

8 

1428 

^(CH^) 

18 

462 

6 (0=0-0) 

9 

1382 

SJOE^) 

19 

406 

6(G-C=0) 

10 

1370 

63(^3) 

20 

299 

6(G-0-C) 

Out- 

-of-plane 

frequencies (a*’ ) 




21 

3020 


26 

712 

y(CH 2) twist 

22 

1428 

S^(OE^) 

27 

583 

y(CH3C02) 

23 

1134 

6(CH^)rock 

28 

238 

Skeletal def. 

24 

948 

yCCH^^wag 

29 

160 

Skeletal def 

25 

874 

y(ch) 

30 

100 

Torsi on (C— CH^ ) 
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3«3.3 Thermodynamic fimctions 

The ideal gas thermal f-unctions C°j, ,S°, 

-(G — Hq)/T have been calculated using the moments of inertia 


obtained from the molecular parameters given in Figure 11 
and the vibrational frequencies given in Table 10. . The 
internal rotational contributions of the methyl, -group to the 
thermodynamic functions have beeh obtained using Pitzer and 
Gwinn's table (62). For this purpose the potential barrier 
height was obtained from equation [19] 


where n = Symmetry numbep of rotating group 
F = h/8 Tt^ cl^ 
h = Planck’s constant 
c = Velocity light 


[19] 


= Reduced moment of inertia which is obtained by 
following formula 


"■ Ai ^ ^i^ 

where = moment of inertia of the rotating group 
= direction cosine 
= principal moment of inertia 

These molecular parameters are given in Table 11. 

As there is some uncertainty in the torsional frequency 
the internal rotational contributions were calculated by varying 
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TABLE 11 ; 

= 47.230 amii 

= 235.820 amu A 
o 

= 279.921 amu A 


molecular pammeters of 

Aa 

2 

2 


= 2 . 390 amu A^ 

O p 

= 2.294 amu At 


f ^ 


c 



7 IRYL ACETATE 
= 0.861 
= 0.509 
= 0.0 

= 100 cm~^ 

= 622.8 cal mol'^ 
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the torsional frequency hy + 10 cm””^. The calculated potential* 
barrier heights are 526, 623 and 728 cal mol~^ corresponding 
to the torsional frequencies 90, 100 and 110 cm~^, respectively. 
The maximum deviations due to the uncertainty in torsional 
frequency are at lower temperatures and they are + 0.08, 

+0.04, +0.05 and +0.09 cal deg“^ mol“^ for C° , , S° 

u 

and -(G-'^-Hq)/T, respectively. This shows that the errors 
in the thermal fimctions are very low owing to the uncertainty 
in the torsional frequency. The enthalpy of formation, G-ibbs 
energy of formation and logarithm of equilibrium constant of 
formation as a function of temperature were Calculated by the 
usual procedure. The enthalpy of formation (63) at 298,15°K 
was taken as -75.46 + 0.18 kcal mol”*^. The values of thermal 
functions of 0(64,65), H2 (66 ), and 82(67) in their reference 
states are used. The calCTilated thermodynamic properties are 
given in Table 04, Appendix 0. 

GaS THERMOSISAMIG PaOPBETISg OF - OROTQHALDEHYPB 
3.4.1 Molecular structure 

Crotonaldehyde can exist in the two forms, S-cis and 
S-trans as shown below. From the microwave (68,69) and 
infrared (70) spectral measurements it has been concluded 
that this molecule- is planar and exists only in the S-trans 
c onf oimat i on , 
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S-cis S-trans 

Su 2 xiki and Kojima (68) and Hsu and Fly gar e (69) have 
independently carried out microwave spectral analysis. They 
reported the moments of inertia and the molecular parameters 
for crotonaldehyde. The agreement between the two sets of 
values is excellent. The molecular structure and structural 
parameters are given in Figure 18. 

3.4.2 Yihrational assignments 

Gredy and Piaux (7l), Saunders et al. (72) and 
Behringer (73) have reported the spectral t®.nds from their 
Raman spectral measurements. Bowles et al. (70) have reported 
almost complete assignments of fundamental frequencies from 
their infrared spectral work in the liquid state. 

Crotonaldehyde belongs to C_ symmetry with nineteen 

o 

in-plane frequencies [19a’] and eight out of plane frequencies 
[Sa'*]. Bowles et al. (70) have reported 25 frequencies. The 
missing ones are two torsional frequencies due to CH^ - and 
CEO - groups. Suzuki and Kojima, and Hsu and Flygare (68,69) 
have reported potential harrier height for CH^- group. Using 
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the equation [19] given earlier in Section 3.3.3, ,the 
torsional frequency due to CH^- group is calculated as 
188 cm*"^. Acraldehyde has similar structure as crotonal- 
dehyde. The torsional frequency due to OHO group is 157 cm“^ 
for acraldehyde. The hand at 300 cm observed in Raman 
spectrum by Saunders et al. (72) for crotonaldehyde could 
be an overtone (2“)^^) of the torsional frequency which 
gives torsional frequency as 150 cm~^. This value has been 
adopted here. The f-undamental frequencies adopted in this 
work are given in Table 12. 

3.4.3 Thermodynamic functions ; 

The ideal gas thermal functions 0°, S°-and 

-(G°-Hq)/T have been calculated using the moments of inertia 
(Table 13) of Suzuki and Kojima (68) and the vibrational 
frequencies given in Table 12, This molecule has two 
rotational groups, CH^- and GH0-. To simplify the cal- 
culations the torsional frequency due to CHO- group is used 
as a fundamental frequency as this would not change the 
thermodynamic values much. The internal rotational contri- 
butions due to CH^-group to the thermodynamic properties 
were obtained using Pitzer and Gwinn's Tables (62). Por 
this purpose, the potential barrier height reported by 
Suzuki and Kojima (68) was used. Reduced moment of inertia (l^) 
was calcula-ted using the parameters given by these authors 
(Table 13) using equation [20], 
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TABLE 12 : ADOPTED FUEDMEM'TAL PREQUM'TCIBS OE OROTOMALDBHYLE 


In-plane frequencies (a’) 


1 

3042 

'^{CH) 

11 

1305 

6(GH) 

2 

3002 

''(CH) 

12 

1253 

6 (OH) 

3 

2944 


13 

1158 

^•(C-G) 

4 

2862 


14 

1109 

6(0H^)rock 

5 

2778 

^Kch(o)) 

15 

1075 

>(C-eH^) 

6 

1693 

V(c=o) 

16 

727 

6( C-CH^) 

7 

1641 

''^(C=C) 

17 

542 

6(GH0) 

8 

1444 

6(CH^) 

18 

459 

6(l:-'C=G-) 

9 

1389 

6(C5H(0)) 

19 

150 

Torsion(C-CHO) 

10 

1375 

6(CH^} 




Out- 

-of-plane 

frequencies (a" ) 




20 

2944 


24 

931 

y(oh(0)) 

21 

1444 

6(01^) 

25 

984 

y(C-GH^) 

22 

1023 

6(CH^)rock 26 

542 

y(::^c=gc) 

23 

966 

Y(aH=GH) 

27 

(188) 

Tor si on (GH^) 

2 
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(Hie enthalpy of formation, A ; the Gibbs energy 
of formation, iiG^ , and logarithm of equilibrinm constant 
of formation as a function of temperature were calculated 
by the usual procedure. The enthalp;^ of formation at 298.15 
was taken as -24.03 kcal mole\rhich is reported by Cox and 
Pilcher (62). The calculated thermodjrnamic properties are 
given in Table C5, Appendix C. 

3.5 AGETOHB. A03TIG ACID AND ACETALDEHYDE 

The thermodynamic functions are available 'in the 
literature for these compounds (66, 74). The enthalpy of 
formation, A ; Gibbs energy of formation, ^ 
logarithm of eqiailibrium constant of formation, log as a 
function of temperature were calcuJ.ated in this work using 
the values of A 298 ^kich are recommended by Cox and 
Pilcher (AhJ 298 -51.9, -103.26, -39.73 kcal mol“^ for 

acetone, acetic acid and acetaldehyde respectively) (63) . The 
thermodynamic properties are given in Tables (T6 - C8, 

Appendix C, 

3.6 ACBTYIME. CARBOM PIOXIDE AFP MTER 

The thermod3mamic properties of ^ 2 ^ 2 ^ ^^2 ^2® 

are available in the literature. The values have been reported 
by TRG (66) and are given in Tables G9 - Gll, Appendix C. 
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3.7 THBRMODTffAMIG PHQPERTIES OP MIX Tim ER 

Tile ■fcliei'niodyiiainic properties of tiie mixture comprising 
^i^iyl acetate, crotonaldehyde, acetone, acetic acid, 
acetaldehyde, acetylene, carbon dioxide and ira.ter are required 
for design of reactor and the separation systems. These 
mixture properties could be calculated for a particular 
composition from the thermodynamic properties of the individual 


components in the mixture (Table S4-C11, appendix G) by the 
following formulae. 

1. Heat capacity of the mixture 


(M) = 




pi 


2. Entha-lpy of the mixture 

- Hg = X X. (H° - Hg) . 

5. Entropy of the mixture 

S? - R X. In X. 
i i ^ X i 

4. Gribbs energy of the mixture 

y , V 

T "4 T li_ + E 1 In 

xi 

where x^ := molar fraction of the individual components. 
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3.8 MTHALPY, G-IBBS ENERGY ASD EQUILIBRIUM CONSTANT OF 
RmCTIOFS 

The following reactions are involved in the production 
of vinyl acetate from acetic acid and acetylene. 

Reactions 

1. CH^COOH + CH^ = CHOOCCH^ 

2. 2CH^C00H — ^ (CE^)^00 + H^O + 00^ 

3- C^B.^ '+ H^O =?- CH^CHO 

4 . 2G2H2 + HgO — GH^OH = GHCHO 

5. n ^ ^*=^2^2^n 

6. CH^COOH + CH^ = GHOOCCH^ s-(CH^G00)2CHCH^ 

The relations for enthalpy, Gibbg energy and equilibrium 
constant, for first four reactions, in terms of temperature 
were derived as follows using the litera,t\ire and the computed 
data (Table G4 - Gll, Appendix C) , These properties were 
calculated for various reaction temperatures and are reported 
in Table G12 - G15, Appendix C. 

3.8.1 Reaction Fumber 1 

GH^GOOH + C 2 H 2 =*► CHg = GHOOCCH^ 

= - 23.457 + 0.02714 T kcal mol”^ 

= - 23.457 kcal mol'‘^ 

InE = lO'* 


- 13.657 
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3*8»2 Reaction Mumber 2 


2GH^C00H 

j 


(CH^)^ CO + 00^ + H^O 


-isg _ _ 96.848 + 8.2 X 10-'^ T + 6.17 x 10'® T%oal nol'l 
= - 96.848 - 6.17 x 10"® teal mol'^ 


_ 48.741 X 10^ , 7 

- T - 4.127 X 10 1 


3.105 X 10^ IT 


3.8.3 Reaction Fumbez* 3 


+ H 2 O 


CH^CHO 


35.951 + 2.382 X lO""^ + 3.750 X 10“^ 

k 'cal niol~^ 

= -35.951 - 3.750 x 10~^ kaal mol“^ 

-i^xr 18.093 X 10^ . 7? 

~ ■ 5 14.254 - 1.887 X 10”^ T 


3 . 8.4 Reaction Number 4 


^^ 2^2 ^ 2 ^ 


CH^CH = CHCHO 


aOj^ - -'75.635 + 6.563 x lo"*^ 1 kical mol~^ 
*^R “ -75.635 kcal 
InK = 10 ^ _ 


+ + + + + 
+ + + 



OOFCLUSION AKD RSCOMMBi® ATI OK'S 
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GOICLUSIOI? 

This study has provided a compre'hensive picture of 
the zinc acetate species on impregnated carbon samples 
both on the surface and in the bulk, as a function of the 
zinc acetate concentration. Evidence is presented for a 
monolayer deposition up to about 15 per cent zinc acetate. 

The zinc acetate snecies, present in the dried catalyst did 
not interact with carrier to form any stable carbon-zinc 
acetate complex at the surface. 

The experimental data for the heterogeneous inter- 
action of acetylene and acetic acid with activated carbon 
based catalysts containing up to 25 per cent zinc acetate 
have been obtained. It is considered that the greatest value 
of this part of the investigation was in setting up a mechanism 
and a rate equation which satisfy the behavior of the acetylene 
-acetic acid system with ZnCOAc)^-^ as. a catalyst. 

The effect of several variables, molar feed ratio of 
acetylene to acetic acid, reaction temperature^ and the reciprocal 
of space velocity on the conversion and product distribution 
were examined. Though fifteen different mechanisms were 
postulated, the i*ate equation, which satisfactorily correlated 
the experimental data, has been evaluated based on the 
mechanism of surface reaction between ionized acetic acid and 
acetylene. The mechanism assumed that the rate controlling 
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step was the irreversible ionic adsorption of acetylene and 
acetic acid. The rate was derived as 


rate(r) 


K 


■A Pa 


1 + 




)(S^) 

Pb 


[ 12 ] 


where and are adsorption equilibrium constants for 
acetic acid and acetylene, respectively. 

Complete vibrational assignments and oalculations of 
the ideal gas thermodynamic properties of vinyl acetate and 
crotonaldehyde hs.ve been made for the first time from 298.15 
to 1500°K using the literature informs-tion on the fundamental 
frequencies and molecular parameters. Internal rotational 
contributions due to CH^-group have been obtained for both the 
molecules using Pitzer and Gwinn's Tables. 

Por vinyl acetate the departure functions, that is, 
the difference between the real gas and the ideal gas property 
have been computed using the Martin-Hou equation of state. 

As the pressure-voliuae-temperature data are not available for 
this compound, the Martin-Hou constants have been predicted 
from the critical constants and the vapor pressure-boiling 
point data of vinyl acetate. These departure functions in 
combination with the ideal gas properties will give the values 
of the real gas thermodynamic properties of vinyl acetate. It 
is hoped that in the absence of the experimental values, these 
predicted values will be of use to the design engineer. 
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REG OHamA-JIOU s 

In the case of reaction studied in this work, additional 
experimental work is needed for further examination of the 
charge transfer (ionic) mechanism. Independent adsorption 
studies of the rea.ction products at reaction conditions along 
with simultaneous measurement of the electron concentration 
may help to determine the electron exchange hetween the 
adsorbed species and the catalyst. 

Experimental pressure-volume-teraperature studies of 
vinyl acetate and the vapor liquid equilibria of the system, 
vinyl acetate-acetic acid-acetene-acetaldehyde_crotonaldehyde 
should be carried out before detailed process engineering 
studies are made. 
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TABLE B12 ; COMP/IRISON OF OBSERVED MB BREBICTED 

MTES AT 155°C 

(¥/E) =30.0 


SI. 

Ex-oerimental 


Predicted 

No. 

PCH^COOH' 

^^2^2 ^ 

Hate, 

"mol gnicat”^ 
iir~^ X 10^ 

Sate , mol 

—1 —1 
gm cat hr 

4 

z 10^ 

1 

0.498 

0..499 

2.10 

2.22 

2 

0.248 

0.749 

1.53 

1.60 

3 

0.198 

0.799 

1.31 

1.35 

4 

0.166 

0.832 

1.15 

1.18 

5 

0.123 

0.874 

0.90 

0.92 


6 


0,099 


0.899 


0.73 


0.75 


112 


SI. 

No. 


1 

2 

3 

4 

5 


M2M_B13: GOMP.iRISON OF OBSERVED aNS PREDICTED 

HITES AT 18Q°G 


( V / E ) = 30.0 


E xperimental 

PCH^COOH’ 


0.496 

0.498 

0.245 

0.749 

0.195 

0.799 

0.163 

0.833 

0.121 

0.874 

0.966 

0.899 


Predicted 

Rate T ^ 

mol gjj.cat~^ Rate, mol 

tix~i X lO'^ ^ cat~-^ lir~^ 
X lo"^ 


4.60 

4.70 

3.90 

4.18 

3.60 

3.69 

3.40 

3.30 

2.93 

2.70 

2.50 

2.25 


6 



■TABLE B14i COMPA: 



SI. Experimental 

t) 2 2 


1 0.494 0.495 

2 0.245 0.748 
5 0.194 0.798 
4 0.162 0.831 
5 0.120 0.874 


6 


0,096 


0.899 
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ION OP OBSERVED ;j?D PREDICTED 


ATES AT 192°G 



Predicted 
Rate, mol 
gm cat~^ lir“^ 

X 10" 


7.12 

6.4 

5.62 

6.1 

4.73 

5.5 

4.10 

4.9 

3.21 

4.1 

2.57 

3.5 


Rate , 

•mol gmcat*"^ 
h.r“^ X lo"^ 





APPENDIX C 


THSEMODYNiiMIC PROPERTIES 



-Ii4 


TABLE CIS ¥APOR PRBSSURS-BOILlK*& POUT DATA FOR ¥IRYL ACETATE 


Antoine equations log mm = 7.2550 - 


1320.2716 


Enthalpy of vaporiC = 8.32 + 0,3 kcal mol' 

zation 


'emperature (t) 


Pressure (p) 
cm Hg 


3.65 

7.80 

13.15 

19.15 
23.05 
27.30 

31.62 
36.94 
42.17 
48.42 

55.62 

61.32 

66. 62 
68.90 
72.10 


3.84 

"4,83 

6.41 

8.68 

10.49 

12.80 

15.59 

19.68 

24.52 

31.56 

41.63 
51.33 
61.91 
66.97 

74.64 
75.63 


72.50 
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[TABLE 


PHYSICAL PROPERTIES OF VINTL ACETATE 


Boiling point, ®C, 760 mm 

72.7 

Freezing point, °C, 760 mm 

-92,8 

Critical temperature, °C 

228.9 

Critical pressure, psia 

609 

Critical volume, litre mol"^ 

0.265 

Flash point, tag open cup,°C 

-5 

Autoignition temperature, ^C 

427 

Explosive limit, lower volume 
per cent in 3,ir 

(M 

Explosive limit, upper volume 
per cent in air 

13.4 


Specific gravity, 20/20°C 0.9338 

Coefficient of cubical 
eznansion ner degree C 

(20-40OC) ■ 1.52 X 10“^ 

Specific heat of liquid, 20°C, 
cal deg~^ 

Refractive index, 

Surface tension, 25°C, dyne 
cja~^ 

Viscosity, 20°C, cp 

Dielectric constant 25°C 

Dipole moment , Debye units 

Themal conductivity at t*^C 
cal sec~^ cm“^ deg~^ 


0.47 

1.3940 

23.8 

0.41 

5.8 

1.75 

370-1.035 (t-20) 
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Table C2 (contd) 



Temperature , 


298.16 

428.16 

438.16 

453.16 

465.16 

472.16 

493.16 


B 

-195. e 

- 438.9 

- 425.3 

- 406.2 

- 392.0 

- 384.0 

- 361.6 



117 


TABLE G3 ; MARTIH-HOU CONSTANTS AND THE DEPARTURE EJNCTIONS 
EOR Vim, ACETATE 




= 0.0 


= 0.082056 


= 0.0, 


Ag 

= -45.20401 

®2 

= 0.040952 


= 0.12129 



= 82.80420 

®3 

= -0.15889 

°3 

=-3851.47 


■^4 

= 2.45173 

=4 

= 0.0 

=4 

= 0.0 


^^5 

= -6.03472 

®5 

=-0.011323 

°5 

= 44.77613 



b = O.i 

33937 

k = -5.0 



T,°E 


P,atm 

7, 



(S-S°), . 




lit. mol cal , mol 


cal mol deg ' 

400 


Saturation -ore ssure = 4,634 atm 





1.60 

28.67 

43.20 


-0.237 



2.31 

21.50 

56.91 


-0.230 



3.09 

17.20 

70.27 


-0.222 



3.94 

14.35 

83.30 


-0.213 

450 


Saturation pressure = 13. 262 atm 





1.55 

28.67 

27.78 


-0.256 



2.16 

21.50 

36.57 


-0.255 



2.80 

17.20 

45.11 


-0.252 



3.46 

14.33 

53.43 


-0.248 



4.14 

12.29 

61.50 


-0.244 



4.83 

10.75 

69.35 


-0.239 



5.51 

9.56 

77.00 


-0.254 



6.19 

8.60 

84.36 


-0.228 



6.85 

7.82 

91.53 


-0.223 



7.50 

7.17 

98.47 


-0.217 
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Table 03 (contd) 


T,°K 


f.atm 

11 aol cal mol"^ 

Saturation pressure = 30.138 atm 




1.57 

28.67 

17.53 

-0.268 

2.14 

21.50 

23.02 

-0.270 

2.71 

17.20 

28.34 

-0.270 

3.28 

14.33 

CO 

-0.270 

3.84 

12.29 

38.45 

-0.268 

4.38 

10.75 

43.26 

-0.267 

4.90 

9.56 

47.90 

-0.264 

5.39 

8.60 

52.37 

-0.262 

5.84 

7.82 

56.69 

-0.259 

6.24 

7.17 

60.85 

-0.256 
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TABLE Oils IDEAL GAS THERMODYNAMIC PRO PERTIES OF WATER 




CONSTANT Af YARIOUS TBIPERATURES 


Reaction No.l 


T,°K 

- A 5° 

k'.oa.l mol~^ 

-^4 

kical mol~^ 

K 

298.15 

15 . 365 

23.457 

1.831 X 10^^ 

428.15 

11.837 

23.457 

1.103 X 10^ 

438.15 

U.567 

23.457 

5.878 X 10^ 

453.15 

11.15$ 

23.457 

2.409 X 10^ 

465.15 

10.833 

23.457 

1.230 X 10^ 

472.15 

10.643 

23.457 

8.447 X 10^ 

493.15 

10.073 

23.457 

2.914 X 10^ 
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TA31B 013 ; EKTMIPY. GIBBS MEIl&Y. AHD ElUI LIBRIUM 
CONSIABT AT V.;lRIOUS TEI'IPERiLTURES 

Reaction No. 2 


T,°K 

.1 

fc.cal mol~ 

- 

iX « 

kcal mol” 

E 

298.15 

96.055 

97.396 

2.529 X 10^^ 

428.15 

95.366 

97. 980 

4.732 X 10^® 

438.15 

95.304 

98.032 

3.412 X 10^'^ 

453.15 

95.209 

98.115 

8.204 X lO"^^ 

465.15 

95.132 

98.183 

4.932 X lO"^"^ 

472.15 

95.085 

98.223 

1.021 X lO'^^ 

493.15 

94.943 

98.348 

42 

1.180 X 10^ 


CONSTilNT AT VARIOUS TSICPERilTURES 


Reaction No. 3 


T°K 

- A 

Real mol”^ 

-A H° ■ 

Real mol”^ 

E 

298.15 

27.174 

36.280 

8.260 X 10^^ 

428.15 

23.137 

36.640 

6.457 X 10^- 

438.15 

22.822 

36.670 

2.399 X 10^^ 

453.15 

22.347 

36.720 

5.984 X 10^® 

465.15 

21.965 

36.760 

2.089 X 10^° 

472.15 

21.743 

36.790 

1.148 X 10^*^ 


2.178 X 10^ 


493.15 


21.072 


56.860 





TABLE 015 ; ENTILILPY. &IBBS. ENERGY, AMB EQUILIBRIUM 
CONSTAMT AT VARIOUS TEMPEELiTURES 

Reaction Ho. 4 


T,°E 

^cal mol 

Real mol 

K 

298.15 

56.067 

. 75.635 

1.230 X lO'^^ 

428.15 

47.535 

75.635 

1.820 X 10^"^ 

438.15 

46.879 

75.635 

2.399 X 10^^ 

453.15 

45.895 

75.635 

1.349 X 10^^ 

465.15 

45.107 

75.635 

1.549 X 10^^ 

472.15 

44.648 

75.635 

4.571 X 10^^ 


19 


493.15 


43.270 


75.635 


1.479 X 10 




131 


COMPUTER PROGlfiJ^ GT 


C 


M4MIM0U COHSTAMS 

i,xx(io) 

SOT:^™%ti6o):v(Jw!if2o6f’®^’®'^’®5'°l’°2,03,C4,C5,B 

COffl'-ION/POINT/lT » V I uo; 

COMOU/liAECI/AK, ill , YS 
COMIyOM/CRITl/TC , PC , 7G ,R , HP 
C OMMOU/POSMi>.T/NIUT , MOUT 

location) .NODKrOR WRITS 

DATaI NIRT,Wt/ 5,6/ 

Mil so, El, E2 /0.2973223s 01,-0.192723845-01,0. 35306520E-04/ 

DATA aIK, AN /5 475, 1.7995/ 

RBxlD (NINT,1)NG0MP 
SORIUT (12) 

NTIME FOR CONTROLLED CARD 
NTIMB=1 

200 WRITE(N0UT,11) 

’ll EORTUT (IHI) 

NK (TITLE(I),I=1,NK) 

SORMiiT (I2,I5A5) 

REilD(NINT,3)TC,PC,7C,B?,R 
PORI-LiT (4F10,4,E12.65 
TB=BP 
GO TO 123 
RKID (NINT,32)N 
S0RMiiT(l3) 

RE.iD(NINT,4)(P(l),I=l,N) 

RK1D(NINT,4)(V(I),I=1,N) 

RE.1D(NINT,4)(T(I),I=1,n) 

PORIC.T (BF10.3) 

70LUI4E COWERSION FROM CC TO LITER 
DO 300 1=1 , N 
7(I)=7(I)/ioOO. 

CONTINUE 
R=R/1000. 

7C=7C/1000. 

YSL=E0+E1‘'TC+E2* TC**2 
YS=EXP ( 2 . 3026 *YSL) 

WRITING OUTPUT FOR INPUT DaIT.! 

WRITE (N0UT,6) (TITLE(I) ,I=1,NK) 

FORJLaT (4X,1545/) 

WRITE (N0UT,7j 

FORMAT (3X,* CALCUIuiTING PLIRTIN-HOU CONSTxlNT BY PREDICTION*) 


2 

3 


122 

32 


300 

123 


6 

7 



o o o 
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W3ITB(I0UT,8) TO PC VO p 

”l?£Uisw>;'s"5:j“ ”/! ’"•”' "iw™ 

» sssr’ **»S. 

GO TO 200 
150 STOP 
MD 

IIBPTC GOYIL 

ColfflOH/MTIN Al ^ ^ 2 ' ^3 * M ^ 

?SS;OT/MmA-'"°°> >'^(200) 

106 |0,PC rO.BP.E 

KPI=-,IB.! ™ RVP10.2,3P8.3P9 4) 

TB=BP 

TBa=BP/TG 

)*--I0&10(P0)-l. 

20=0. 291-0. 080*¥ 

C^iLCUL.lTIFG AG 

TBofM.*S?4°*Toff 6?5 |o ™pf^fraB,DBG.K: 

CTn?j^ *!?• 9869-0. 6751*ZC ) 

AfiJro^sf “>‘^3 B 

DBPIBB V-iLGE — ill Yq 

iiMS=-.m*PC/TG ’ 

B=il *!rc *( ZC-BTA /15 ) /PG 
iiN=1.70 /-^:5.;/±'o 

YS=3.0 

M/iPTIN-HOU PAiLAiETERS 

31=R 

ci=0 

K=12. *P0*(vis)24-3:"A*^?0* 
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C25=( (Y2+B*R*!rpi + (H*-TFI)**2/PC) *(’1 

.?l:ii3;."!;j«Hig5‘"'"™''”|-"‘'^“»)'i™-’:i 

(f^~^)f*^-((yc/AiN)-B) *^^)/( (VC/UJ-31^f*2-fvn 

(< -(vaA.-3)**2) 

B3=B£i— BxLt,.-.BiLiA 

ii3=73_B3*TG~C3*SX;p (KFI) 

i^4=T4 

B4=0. 

C4=0. 

C|=^C2*(7C-B)**5-C.3je (TC-B)**2 

«3-B3^TC-B)«2 

WRIJE OUTPUT 

form! T ^ ( lOX ! 2H?it^TRY ’ «l ^ 1 » 1^! , B5 , C5 



100 


32HC5/3B20,8) 

WRITE (ROUT, 105)B,KEI 
RSTOffi ^^®’1™>1SX,2HK/3B20.S/) 
BED 


lEETRT 



oooaaooaooooooa aaoo 
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COMPUTER PRQgBAM 02 


THIS PROGRAM USES MARTIH-HOU EQUATION AND MARTIN 
EQUATION (T.P.)TO CALCULATE THERjyiAL FUNCTION AND P¥T 
DATA 

MARTIN-HOU EQUATION IS P=RT/(V-B)+SIGMA(A(I)+B(I)*T+C(I) 
1*EXP(AK*T/TC) )/(¥-B) 

P(PSI),T(DBG.R),¥ AND B(FT3/LB). 

A2..A5,B2..B5,C2..C5 ARE CONSTANT IN MARTIN-HOU EQUATION 
R=10. 7351, (GAS CONSTANT) , TO (DEG. R) . 

TEMPI, TMIN,TMAX,TLIMIT,TRIS3 ALL IN DBG R 
MARTIN EQUATION IS . .L0GP(ATM)=X1+X2/T+X3*T+X4*T 2+X5* 

T 3+L0G10(TC-T+8.0) T IN DEG.E. 

DM IS DENSITY, TEMP IS TSI'-IPBaATURE IN DEG K 
THIS PROGRAM CALCULATES THERMAL FUITCTIOHS IF DENSITY 
IS KNOWN AT ANY TMPERATURE 

TMX=UPPER LIMIT FOR TMP. . , TMIN=LOWER LIMIT FOR TMP. 
TRISE=INCREMMT AT WHICH P¥T Y-iLUES ARB WANTED 
TMPI=INITIAL TE€P., ¥M=M01.¥T. 

SEQUMGS OF DATA CARDS. .TBI''IPI,TMIN,TI4AX,TLIMIT,TRISS 
{6F10.5),A2,B2,C2,...A5,B5,C5,TG,M,R,B,AK{6F10.5). , 
X1,..X6(6F10.6).,TBMP <F10.5). 

DIMMSION DEN (60), PGALC(60) 

PRINT 50 
PRINT 20 

READ 9,TEMPI,TMIN,TMAX.TLIMrT,TRISE 
READ 8,A2,B2,C2. 

READ 8,A3,B3,C3, 

READ 8,A4,B4,C4 
READ 8,A5,B5,C5 

8 FORM.iT (E15.8,4X,E15.8,5X,E15.8) 

REiiD 7, TC,¥M,R,PC 

7 FORMAT (F8.3,F8.3,F12.6,F8.2) 

TC=TC *1.8 

C TC MUST BE FED IN DEGREES KEL¥IN 

R=R/WM 
TEMP=TEMPI 

9 FORIAIT (6F10.5) 

10 FORNIiiT (F10.5) 

DEN.(L)==0,Q03 

N=10 

DO 11 1=1, N 
11 = 1+1 

DENt 1 I ) =DBN (I ) +0 . 001 

11 CONTINUE 
AK=5.0/TC 

100 PI=0.0001 

PRINT 40 
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IJCTBI'iP.GE. TMI l. AND, 'TEMP. LT.m^^)RBiiD 10,TMP 
T=TEMP 

IP(T*1.8.GT.TG) GO TO 101 

C THE ¥.iLaE OP SHOULD BE IN ..TMOSPHERBS 

PMM=7. 255-1320. 2716/(T-45. 96) 

PP=(10.0 *«»MM)/760.0 
PHIZTT 1,PP 

FORiUT(/55X,* aiTUiLiTIOH PEBSSiniE=*, S15.6) 

101 EX=HXP (-AK*T*I . 8 ) 

PCALC(1)=0.0 

PORTILIH SOURCE LIST FUIN 

ISN SOURCE STATEMENT 

PC=PC*14.696 
B=R«C/(8.0*PC) 

T=TMP*1.8 

DO 60 1=1, M 

V=1.0/(DBN(I)*62.482) 

VB=V-B 

P1T=H*T/¥B 

P2T= ( A2+B2 *T+C2*EK ) /¥B «t2 
F5T= ( A3+B3 ^+03 *EX ) /VB **3 
F4T= ( A4+B4 *T+C4*EX ) /¥B**4 
P5T=(A5+B5*T+C5*BX)/¥B**5 
PG..LC ( I+l )=FlT-rF2T+F3T+F4T+F5T 
IF(T.GE.TG) GO TO 14 
IP(PC..LC(I+1).&T.PP*14.7) GO TO 60 
IF (aBS (PGxiLC(l+l)),LT.ABS(PC..LG(l))) GO TO 88 
14 P=PG-iLC(I+l) 

TERM1=P*¥-R*T 

TERM2 = (A2+{1.0+iIK;*T)*C2*EX)/¥B 

TERM3 = (A3+(1.0+.IK*T)*C3*1X)/{2.0 ¥B**2) 

TERM4 = (il4+(l.0+iiE:*T)*G4*BX)/(5.0 ¥B»*3) 

TERM5 = (A5+(l.0+AK*T)*C5*EX)/(4.0 ¥B**4) 
DBLH=TERMl+TERM2+TBRM3+TERM4+TERM5 
¥BB=R*T/P(I)-B 
TERMl=R5f ALOG ( ¥B/¥EB ) 

TBRM2=- ( B2-aE^ G2*BK ) * ( 1 . 0/¥B-l . 0 ABB ) 

TERM3=-0 . 5* ( B3-AK*C3 *EX ) *( 1 . 0/¥B**2-l . 0 ABB»*2 ) 

TERM4=-1 . 0/3 . 0* ( B4-.'1E*04 *EX ) * ( 1 . 0/¥B** 3-1 . 0/¥BB **3 ) 

TERM5=-0.25*(B5-.iE:C5*iS;)*(l.O/¥B#*4-1.0/¥BB *4) 

DELS=TBRM1+TERM2+TERM3+TERM4+TBRI!5 

DELH=DELH*1 . 98589/10. 7315 

DELS=DELS*1 . 98589/10 . 7315 

TBMPK=T/1.8 

PATM=P/I4.f96 

¥CGS=VM/DEN(I) 

PRINT 30, TMPZ, PATH, ¥CGS, DELS, DELS 
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60 CONTII'Ttjb 
88 COKTIRUS 

IF ( TS4P , OE. TMIN.iiND. IMP . LT. miK ) GOTO 100 
TMP=TEMP+TRISE 
IF (TEKP.BQ.TLIMIT) STOP 
00 TO 100 


20 FORIvLiT (50X,87(1H*)//30X,87H *TBI©BiLlTURE * pressure 
1* VOLUME * DELTA-«^H * DELTii S * /30X, 8 7H * DEGREES * UTMOS. 
_ 2 -GC2GM MOLE*C.il/GM.MOLB* CAL/(GH.MOLS, DEGREES)*/ ) ' 

4X, F9.3, 4X,1H*, 4X,E15.6,3X,1H*, IX, 

50 FORIlxT Uhi) 

70 COETIITUE 


EED 


^ITTRY 


